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Hydraulic  and  Physical  Properties  Affecting  Ice  Jams 

KATHLEEN  D.  WHITE 


INTRODUCTION 

River  ice  jams  can  form  suddenly,  bringing 
about  rapid  fluctuations  in  stage.  Ice  jams  can 
cause  flooding  upstream  and  decrease  down¬ 
stream  discharge.  Damage  to  riverine  structures 
such  as  bridges,  locks,  dams,  dikes,  groins,  levees, 
and  riprap  installations  is  also  possible.  Ice  jams 
affect  navigation  through  delays,  stoppages,  or 
damage  to  tows,  barges,  and  mooring/fleeting 
areas.  Ice  can  block  hydropower  and  water  sup¬ 
ply  intakes.  Ice-induced  scour  may  cause  the  ero¬ 
sion  of  river  beds  and  banks,  with  adverse  impacts 
on  fish  and  wildlife  habitat,  as  well  as  the  expo¬ 
sure  of  utilities  buried  beneath  the  streambed. 
Emergency  and  medical  relief  to  flooded  areas 
may  be  limited  by  flooding  or  ice-related  scour  and 
erosion  of  roads  resulting  in  road  closures,  or  by 
the  closure  of  bridges  weakened  or  destroyed  by 
ice.  The  potential  exists  for  death  or  serious  injury 
during  ice-related  flooding,  evacuations,  and  other 
ice  mitigation  operations. 

Ice-related  damages  can  be  minimized  or 
avoided  through  the  use  of  hydraulic  modeling 
to  predict  jam  location,  upstream  and  downstream 
stages,  jam  thickness,  and  other  jam  characteris¬ 
tics.  While  knowledge  about  ice  jam  formation, 
progression,  and  failure  has  increased  in  recent 
years,  to  date  no  deterministic  model  of  these  ice 
jam  processes  has  been  developed.  In  part,  these 
complex  physical  processes  are  poorly  understood 
because  field  observations  of  ice  processes  and 
properties  are  spatially  and  temporally  limited  and 
can  be  dangerous  to  acquire.  Data  are  often  scat¬ 
tered,  and  although  earlier  literature  reviews  (e.g., 
Bolsenga  1968)  contain  a  great  deal  of  historical 
information  regarding  river  ice  covers,  they  do  not 
include  the  tremendous  amount  of  information 
gathered  over  the  past  few  decades.  As  a  result, 
many  of  the  variables  and  parameters  included  in 
the  empirical,  analytical,  and  numerical  models 


that  do  exist  for  some  ice  jam  processes  are  not 
sufficiently  well  known  or  described. 

Beltaos  (1995)  and  Ashton  (1986)  provide  thor¬ 
ough  discussions  of  the  various  processes  impor¬ 
tant  in  ice  cover  initiation,  growth,  breakup,  trans¬ 
port,  and  jamming.  The  purpose  of  the  present 
study  is  to  review  available  information  on  the 
properties  of  ice  covers  and  ice  jams  that  are  im¬ 
portant  in  modeling  river  ice  jam  hydraulics.  An 
overview  of  ice  jam  hydraulics  identifies  specific 
ice  properties  of  interest,  which  are  discussed  in 
detail  in  the  following  section.  Information  is  also 
presented  about  miscellaneous  ice  properties  ad¬ 
dressed  in  the  literature  but  not  explicitly  included 
in  hydraulic  models  of  ice  jams. 

OVERVIEW  OF  ICE  JAM  HYDRAULICS 

Formation  processes  of  river  ice  covers  range 
from  the  purely  static  to  the  purely  dynamic.  Static 
ice  cover  formation  is  a  largely  thermal  process  in 
that  the  initiation  and  growth  of  ice  crystals  result 
from  heat  transfer  between  the  water  and  the  at¬ 
mosphere.  Statically  initiated  ice  covers  in  rivers 
are  found  in  quiescent  areas  and  along  the  banks 
(border  ice).  Dynamic  ice  cover  formation  results 
from  the  mechanical  processes  associated  with  floe 
interactions.  These  may  range  from  relatively  low- 
energy  processes  such  as  the  juxtaposition  of  ice 
floes  into  a  single  layer  of  ice  that  then  freezes  in 
place,  or  higher-energy  processes  such  as  the  ac¬ 
cumulation  of  floes  into  an  ice  jam  by  shoving  and 
internal  collapse.  Temperature  plays  an  important 
role  in  ice  cover  formation.  Andres  (1999)  states 
that  for  large,  moderately  sloping  rivers,  colder 
temperatures  will  result  in  juxtaposition,  while 
higher  temperatures  will  result  in  ice  covers 
formed  by  shoving  and  internal  collapse  (or  "con¬ 
solidation").  However,  he  finds  that  for  very  mild 
slopes  and  very  steep  slopes,  the  process  is  rela- 


tively  independent  of  temperature:  if  the  river 
slope  is  mild,  juxtaposed  ice  covers  form,  and  if 
the  slope  is  steep,  the  ice  cover  will  form  from 
shoving.  Once  formed,  ice  covers  can  also  thicken 
via  thermal  processes,  by  flooding  and  refreezing 
of  the  surface,  or  by  transport  and  deposition  of 
ice  beneath  the  cover.  An  overview  of  the  processes 
resulting  in  the  formation  of  ice  covers  and  ice  jams 
is  presented  below. 

Ice  cover  formed  from  juxtaposition  of  ice  floes 
Juxtaposition  occurs  when  frazil  or  solid  ice 
floes  carried  downstream  reach  an  obstacle  (such 
as  an  ice  cover)  and  stop,  forming  a  single  layer 
ice  cover  (Fig.  1).  Freezeup  accumulation  and  pro¬ 
gression  due  to  juxtaposition  can  be  quite  rapid. 
For  example,  Gerard  and  Andres  (1982)  noted  rates 
of  up  to  7  km/day  for  a  freezeup  jam  on  the  North 
Saskatchewan  River,  Alberta,  while  Beltaos  (1983) 
reported  8.4  km/day  over  a  20-hour  period  on  the 
Thames  River,  New  Brunswick.  The  physical  pro¬ 
cess  of  the  juxtaposition  of  frazil  or  solid  ice  floes 
into  a  single  layer  ice  cover  has  not  yet  been  care¬ 
fully  described,  perhaps  because  these  ice  accu¬ 
mulations  are  less  threatening  than  thickened 
jams.  Yet,  juxtaposition  is  often  the  primary  pro¬ 
gression  mode  for  initial  ice  cover  formation  in 
large  rivers  and  for  freezeup  jams.  Traditionally, 
two  approaches  have  been  taken  in  determining 
whether  the  incoming  ice  floes  will  be  stable  and 
juxtaposition  is  possible:  use  of  a  surface  velocity 
criterion  or  of  a  criterion  describing  the  rotational 
stability  of  the  incoming  ice  floes. 


The  critical  velocity  would  be  the  maximum  ve¬ 
locity  at  which  juxtaposition  can  still  occur;  higher 
velocities  would  result  in  underturning  of  ice  floes 
or  shoving.  Majewski  and  Grzes  (1986),  for  ex¬ 
ample,  reported  that  velocities  less  than  about  0.5 
to  0.7  m/ s  allowed  juxtaposition  in  one  reach  of 
the  Vistula  River,  Poland,  while  greater  than  0.6 
to  0.7  m/s  caused  shoving.  The  idea  of  a  critical 
velocity  criterion  is  attractive,  yet  it  fails  to  account 
for  many  of  the  important  variables,  such  as  ice 
discharge,  and  the  thickness,  size,  and  cohesion 
of  the  ice  floe.  As  a  result,  a  surface  velocity  crite¬ 
rion  is  not  acceptable  except  in  limited,  site- 
specific  cases. 

Efforts  to  expand  the  critical  surface  criterion 
led  to  variations  on  the  Froude  number  (e.g., 
Pariset  and  Hausser  1961,  Pariset  et  al.  1966).  The 
Froude  number,  F,  is  the  ratio  of  inertial  forces  to 
gravity  forces: 

&  m 

where  v  is  the  average  velocity  upstream  from  the 
ice  cover,  g  is  gravity  (981  cm/s2),  and  d  is  some 
appropriate  length  scale,  generally  taken  to  be  the 
upstream  flow  depth,  H,  for  open-water  situations. 
A  more  appropriate  length  scale  when  consider¬ 
ing  ice  transport  would  be  the  ice  piece  thickness. 
For  ease  in  application,  a  juxtaposition  criterion 
based  on  F  is  often  presented  in  the  form  of  a  ve¬ 
locity  criterion.  Pariset  and  Hausser  (1961)  pre¬ 
sented  their  Froude  criterion  in  the  form  of  a  criti¬ 
cal  velocity,  but  as  they  observed,  ice  covers  have 


Figure  1.  Typical  juxtaposed  ice  cover ,  formed  here  from  large 
broken  pieces  of  sheet  ice . 
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been  known  to  progress  by  juxtaposition  at  veloci¬ 
ties  more  than  three  times  greater  than  the  critical 
velocity  they  predicted.  Using  a  different  ap¬ 
proach,  Tremblay  and  Thibeault  (1980)  proposed 
a  critical  velocity  criterion  for  the  Milles  lies  River 
(Quebec)  based  on  the  Chezy  formula  that  in¬ 
cluded  floe  size  and  ice  discharge. 

Later  researchers  (e.g.,  Uzuner  and  Kennedy 
1972, 1974;  Ashton  1974;  Uzuner  1975;  Michel  1978; 
Tatinclaux  and  Gogus  1981)  used  both  laboratory 
experiments  and  analysis  to  arrive  at  a  rational 
Froude  number  to  describe  the  limits  of  the  rota¬ 
tional  stability  of  the  ice  floes.  Generally,  a  densi- 
metric  Froude  number  was  used,  with  a  length 
scale  based  on  some  ratio  between  the  ice  block 
thickness  and  water  depth,  modified  by  the  dif¬ 
ference  in  the  densities  of  water,  p,  and  ice,  p*. 
Ashton  (1974)  successfully  described  the  experi¬ 
mental  data  using  an  analysis  based  on  the  "no- 
spill"  condition,  which  postulates  that  the  limit  of 
rotational  stability  occurs  when  the  upper,  up¬ 
stream  comer  of  the  block  becomes  submerged. 
He  assumed  that  the  underturning  moment  was 
caused  by  the  acceleration  of  flow  under  the  block 
because  of  the  constriction  caused  by  the  block. 
Despite  the  obvious  problem  that  the  constriction 
of  flow  would  be  negligible  in  deep  water, 
Ashton's  criterion  was  widely  adopted. 

Daly  and  Axelson  (1990)  demonstrated  that 
floating  blocks  have  a  well-defined  maximum 
righting  moment.  They  noted  that  the  limit  of  ro¬ 
tational  stability  has  been  reached  when  the 
underturning  moments  exceed  this  maximum  and 
a  no-spill  condition  was  not  necessary.  They  de¬ 
veloped  an  analytical  expression  for  the  limit  of 
rotational  stability  in  terms  of  the  hydrostatic  right¬ 
ing  moment  of  the  block,  expressed  in  the  form  of 
a  densimetric  Froude  number.  Their  work  was 
incorporated  by  McGilvary  and  Coutermarsh 
(1992)  in  a  dynamic  analysis  of  ice  block  stability 
based  on  careful  measurements  of  the  pressure 
distributions  around  the  leading  edge  and  under¬ 
side  of  a  model  ice  floe  for  various  flow  depths 
(Coutermarsh  and  McGilvary  1991,  1993).  This 
pressure  distribution,  produced  by  the  accelera¬ 
tion  of  the  flow  as  it  passed  under  the  block,  cre¬ 
ates  an  underturning  moment  acting  on  the  block 
that  could  cause  undertuming  in  both  deep  and 
shallow  water. 

McGilvary  and  Coutermarsh  (1992)  further  ex¬ 
tended  the  analysis  by  including  the  angular  mo¬ 
mentum  of  the  block,  and  demonstrated  that  the 
angular  momentum  could  be  an  important  factor 
in  causing  blocks  to  underturn  (all  of  the  previous 


analyses  had  been  static  in  that  they  neglected  the 
angular  momentum  of  the  rotating  blocks).  They 
found  that  the  Froude  criteria  of  Pariset  and 
Hausser  (1961)  and  Ashton  (1974)  underpredicted 
the  critical  Froude  number  for  stability,  particu¬ 
larly  for  small  ratios  of  floe  thickness  to  depth. 
They  did  find  good  agreement,  however,  with  the 
shallow  water  criterion  for  stability  presented  by 
Daly  (1984): 


where  t ^  is  ice  block  thickness,  L  is  ice  block  length, 
H  is  upstream  water  depth,  and  P  is  a  coefficient 
that  is  independent  of  underturning  velocity,  and 
block  length,  thickness,  and  density,  determined 
to  be  equal  to  0.262.  A  definition  sketch  is  given  in 
Figure  2.  The  submerged  depth  of  the  ice  block,  ts, 
is  calculated  using 


t 


s  — 


(3) 


Daly  (1984)  noted  that  this  equation,  developed 
as  a  simple  curve  fit  to  data  from  Uzuner  and 
Kennedy  (1972),  Ashton  (1974),  and  Larsen  (1975), 
is  not  strictly  applicable  to  deep  water  (ice  block 
thickness-to-depth  ratios  less  than  about  0.1). 

More  recently,  Andres  (1999)  has  proposed  the 
use  of  a  dimensionless  stability  number  to  deter¬ 
mine  whether  juxtaposition  or  shoving  will  domi¬ 
nate  the  ice  cover  formation  process  under  differ¬ 
ent  hydrologic  and  meteorological  conditions.  His 
analysis  includes  the  effects  of  temperature  in  the 
form  of  air  temperature  (Ta),  the  thermal  conduc¬ 
tivity  of  the  ice  (/q),  and  the  latent  heat  of  fusion  of 
the  water  (Lf).  The  dimensionless  stability  num¬ 
ber  takes  the  form 


X  = 

QSfpiLf 


(4) 


where  Q  is  discharge,  B  is  river  width,  and  Sf  is 
the  water  slope.  A  typical  value  of  the  latent  heat 
of  fusion  of  water  is  334  J/g  (Batchelor  1967).  The 
thermal  conductivity  of  the  ice  is  related  to  tem¬ 
perature  by  the  following  equation  (US ACE  1999): 
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^=2.21-0.117*  (5) 

where  T*  is  the  temperature  of  the  ice  (°C)  and  /q 
has  units  J/m  s  °C.  Andres  reported  that  the  criti¬ 
cal  dimensionless  stability  number  separating  jux¬ 
taposition  and  shoving  on  the  Peace  River  in  north¬ 
ern  Alberta  is  about  0.0003,  with  shoving  more 
often  observed  at  lower  values  and  juxtaposition 
most  often  observed  at  higher  values.  He  noted 
that  river  slope  and  crushing  strength  of  the  ice 
(again,  related  to  temperature)  contributed  to  the 
uncertainty  in  determining  the  critical  dimension¬ 
less  stability  number.  Although  not  widely  tested 
at  this  time,  this  method  holds  promise  for  deter¬ 
mining  juxtaposition  potential  for  ice  cover  for¬ 
mation  and  freezeup  jams. 

HEC-RAS  assumes  single-layer  juxtaposition  of 
floes  if  the  ice-jam  force  balance  (see  section  on 
shoving  and  internal  collapse  below)  shows  that 
the  thickness  of  the  parent  ice  cover  making  up 
the  floes  exceeds  that  necessary  to  balance  forces 
in  the  downstream  direction.  The  ICETHK  option 
of  the  USACE  step-backwater  computer  program 
HEC-2  (USACE  1990)  uses  the  criteria  developed 
by  Ashton  (1974)  and  Michel  (1978)  to  determine 
when  juxtaposition  may  occur.  However,  it  is  rec¬ 
ommended  here  that  eq  2  be  applied  to  variables 
derived  from  ICETHK  and  other  numerical  mod¬ 
els  for  the  determination  of  juxtaposition  poten¬ 
tial  because  of  its  good  agreement  with  the  results 
of  the  fairly  complex  model  presented  by 
McGilvary  and  Coutermarsh  (1992).  Application 
of  eq  2  will  require  some  knowledge  or  estima¬ 
tion  of  the  following  ice  properties:  t^,  p,  ,  and 
H,  of  which  two,  p  and  p±,  are  well  known  (1000 
kg/m3  and  916  kg/m3,  respectively). 


Ice  cover  formed  by  undercover  deposition 

Undercover  deposition  of  ice  occurs  when  in¬ 
coming  frazil  or  solid  ice  floes  are  unstable  at  the 
upstream  end  of  an  ice  cover  (e.g.,  juxtaposition 
criterion  is  exceeded).  Once  the  floes  underturn, 
they  can  be  transported  beneath  the  ice  cover  and 
into  the  open  water  beyond,  or  they  may  deposit 
beneath  the  ice,  depending  on  the  river  hydrau¬ 
lics  and  the  ice  properties  (Fig.  3).  Predictions  of 
where  and  how  much  under-ice  deposition  will 
occur  have  been  made  using  three  approaches: 
critical  velocity;  Froude-number-based  block  sta¬ 
bility  criteria  similar  to  those  discussed  above,  but 
for  submerged  ice  floes;  and,  more  recently,  ap¬ 
plication  of  bed-load  transport  theory.  Ice  jam 
thickening  by  deposition  of  ice  is  not  specified  in 
either  HEC-RAS  or  ICETHK,  although  ICETHK 
does  allow  thinning  of  an  ice  cover  by  erosion 
when  the  critical  velocity  is  exceeded. 

The  critical  deposition  velocity  approach  as¬ 
sumes  that  frazil  deposits  (or  hanging  dams)  ac¬ 
cumulate  from  upstream  to  downstream,  begin¬ 
ning  with  the  most  upstream  location  where  the 
average  cross-sectional  flow  velocity  is  less  than 
some  critical  deposition  velocity.  Frazil  will  de¬ 
posit  at  that  location  until  the  critical  velocity  has 
been  reached.  Any  additional  frazil  will  be  trans¬ 
ported  downstream  to  some  point  where  the  ve¬ 
locity  is  less  than  the  critical  velocity,  where  the 
process  begins  again.  According  to  Beltaos  (1995), 
this  critical  velocity  method  is  used  in  the  numeri¬ 
cal  model  ICESIM  to  estimate  the  location  and 
thickness  of  frazil  deposition. 

Observations  reported  by  Michel  and  Drouin 
(1981)  indicate  that  critical  velocity  can  vary  over 
the  course  of  a  winter.  They  reported  that  the  value 
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Figure  3.  Frazil  deposit  formation. 


decreased  from  0.9  m/s  in  early  winter  to  0.6  m/s 
in  late  winter  on  the  LaGrande  River,  Quebec.  Daly 
(1994)  attributes  reported  variations  in  critical  ve¬ 
locity  to  factors  such  as  the  form  of  the  frazil,  the 
length  of  travel,  and  the  degree  of  supercooling  in 
the  river.  Michel  (1984)  suggests  that  variations  in 
critical  velocity  may  be  explained  by  the  physical 
state  of  the  frazil  particles.  In  his  experience,  criti¬ 
cal  velocities  decrease  with  age  from  a  high  of 
about  1.4  m/s  for  active,  or  recently  formed,  frazil 
to  a  low  of  about  0.6  m/s  for  aged  frazil  floes. 

These  variations,  along  with  variations  in  the 
critical  velocity  between  rivers  and  between  dif¬ 
ferent  reaches  on  the  same  river,  made  it  difficult 
to  apply  this  method.  In  view  of  these  difficulties. 
Sun  and  Shen  (1988)  proposed  an  empirical 
method  to  estimate  ice  thickness  based  on  Froude 
number,  similar  to  the  approach  taken  for  juxta¬ 
position.  The  stability  criteria  resulting  from  labo¬ 
ratory  studies  and  analytical  analyses  (e.g.,  Ashton 
1974,  Gogus  and  Tatinclaux  1981,  Daly  and 
Axelson  1990)  have  taken  the  form  of  some  varia¬ 
tion  on  the  densimetric  Froude  number.  However, 
Starosolzky  (1981)  suggested  that  scatter  in  experi¬ 
mental  data  for  parallelepipeds  using  this  method 
could  be  explained  by  a  bed-load  transport  ap¬ 
proach. 

Although  not  yet  applied  to  blocks,  bed  load 
transport  has  been  applied  to  frazil  ice  by  Shen 
and  Wang  (1992),  who  proposed  a  new  method  to 
analyze  frazil  transport  based  on  the  concept  that 
the  movement  of  the  ice  particles  is  caused  by  the 
shear  stress  acting  on  the  particles.  Essentially,  the 
method  relies  on  the  functional  relationship  be¬ 
tween  the  dimensionless  ice  transport  capacity,  <t>, 
and  the  dimensionless  flow  strength,  0: 

O  =  /(©).  (6) 


The  dimensionless  ice  transport  capacity  is  de¬ 
scribed  by 


in  which  qi  is  the  ice  transport  capacity  per  unit 
width  and  d  is  the  particle  diameter.  The  dimen¬ 
sionless  flow  strength  is  calculated  using 

0_(p-pi)^  (8) 

where  t  is  the  shear  force  on  the  underside  of 
the  ice  cover.  Laboratory  studies  reported  by  Shen 
and  Wang  (1995)  suggested  that  the  dimension¬ 
less  parameters  and  0  in  eq  6  and  7  should 
be  modified  to  account  for  the  effect  of  the  shape 
of  the  particles,  resulting  in  the  following  formu¬ 
lations: 


where  Cf  is  a  fall  velocity  coefficient,  dn  is  the  nomi¬ 
nal  particle  diameter,  and  u*  is  the  shear  velocity: 


5 


where  R j  is  the  hydraulic  radius  of  the  ice-covered 
channel,  and  t  is  the  shear  stress  on  the  underside 
of  the  ice  cover. 

Because  the  width  of  a  river  channel  is  often 
much  wider  than  its  depth,  the  open-water  hy¬ 
draulic  radius  R  (the  ratio  of  cross-sectional  area 
to  wetted  perimeter)  of  an  open-water  reach  is  of¬ 
ten  approximated  as  H,  the  water  depth.  In  the 
ice-covered  case,  the  wetted  perimeter  is  effec¬ 
tively  doubled,  so  Rj  can  be  approximated  by  y-J 
2,  where  y±  is  the  under-ice  depth,  often  calculated 
from  Manning's  formula  as 


Vi= 


(12) 


where  c  is  a  coefficient  (1.5875  in  SI  units  of  s/m1/3 
and  1.0654  in  English  units  of  s/ft1/3),  q  is  the  dis¬ 
charge  per  unit  width,  and  nc  is  the  composite 
Manning's  n  for  the  channel.  Nezhikhovskiy 
(1964)  reviewed  methods  of  determining  compos¬ 
ite  channel  roughness  and  suggested  the  use  of 
the  Belokon-Sabaneev  formula: 


(13) 


in  which  n f  is  the  Manning's  n  for  underside  of  ice 
and  ttb  is  the  Manning's  n  for  the  bed.  Other  for¬ 
mulations  are  presented  in  Uzuner  (1975).  The 
roughness  ratio  can  be  used  to  modify  the  ice- 
affected  hydraulic  radius  as  follows: 


V\_ 
2  ‘ 


(14) 


White  and  Acone  (1998)  modeled  frazil  depo¬ 
sition  beneath  an  ice  cover  at  a  river-reservoir 
confluence  utilizing  both  ICETHK  and  the  sim¬ 
plified  version  of  Shen  and  Wang's  functional  re¬ 
lationship  in  eq  6  to  10: 


4>n=5.487(0n-0c)1-5  (15) 

where  0C  =  0.041.  The  following  process  was  used: 
1)  a  hydraulic  model  with  stable  ice  cover  must 
be  developed  for  the  reach  of  river  in  question;  2) 
variables  R  and  Sf  calculated  by  the  hydraulic 
model  are  used  to  estimate  0n  using  eq  10  and  11, 
from  which  On  is  calculated  using  eq  15;  3)  frazil 


transport  is  assumed  to  occur  at  locations  where  a 
solution  to  eq  15  is  obtained;  if  no  solution  is  pos¬ 
sible  (i.e.,  On  <  <I>C),  deposition  is  assumed  to  oc¬ 
cur.  An  iterative  process  is  required  to  obtain  an 
equilibrium  frazil  deposit.  The  ice  cover  of  the 
hydraulic  model  is  thickened  at  the  locations 
where  deposition  is  identified,  beginning  with  the 
most  upstream  location.  The  ice-covered  hydrau¬ 
lic  model  is  rerun,  providing  new  values  of  R  and 
Sf,  which  are  used  to  calculate  a  new  On  for  each 
cross  section,  thus  identifying  locations  of  trans¬ 
port  and  deposition  under  the  updated  hydraulic 
conditions.  The  process  is  repeated  until  no  fur¬ 
ther  deposition  is  possible  or  until  the  volume  of 
the  deposited  frazil  ice  nears  the  estimated  vol¬ 
ume  produced.  The  variables  required  in  this 
method  are  p,  pir  dn,  Cf,  Sf,  and  the  variables  nec¬ 
essary  to  calculate  R:  q,  nif  n^,  and  nc. 

Ice  cover  formed  from  shoving  and 
internal  collapse 

When  the  forces  exerted  by  incoming  ice  in¬ 
crease  to  the  point  that  they  exceed  the  internal 
strength  of  ice  accumulations  formed  by  juxtapo¬ 
sition  or  thickened  by  deposition,  the  ice  accumu¬ 
lation  will  fail.  The  failure  may  be  large  scale,  re¬ 
sulting  in  the  transport  of  ice  for  large  distances 
until  the  accumulation  stops  again,  or  it  may  be 
smaller  in  scale,  resulting  in  a  local  shoving  and 
thickening  of  the  accumulation.  The  shoving  and 
thickening  process  is  an  inherently  unsteady  pro¬ 
cess  that  may  eventually  cause  the  formation  of 
what  is  often  termed  an  equilibrium  jam;  that  is,  a 
jam  in  which  a  uniform  section  of  ice  has  reached 
a  thickness  where  the  external  and  internal  forces 
are  in  equilibrium  (Fig.  4). 

Zufelt  and  Ettema  (1997)  developed  an  un¬ 
steady  model  of  the  shoving  and  thickening  pro¬ 
cess,  and  a  stationary  ice  cover  subroutine  has  been 
used  in  conjunction  with  the  one-dimensional 
unsteady  flow  model  UNET  developed  by  the 
Hydrologic  Engineering  Center  (HEC)  (USACE 
1995).  The  use  of  a  dynamic  ice  cover  in  UNET  is 
under  development  (McGilvary  et  al.  1995,  Daly 
et  al.  1997).  However,  most  currently  used  mod¬ 
els  assume  steady  state  conditions,  including  the 
HEC's  River  Analysis  System  (HEC-RAS)  (USACE 
1997),  the  ICETHK  option  of  the  USACE  step- 
backwater  computer  program  HEC-2  (USACE 
1990),  and  Canadian  models  RIVJAM  and  ICESIM 
(Healy  et  al.  1997).  This  assumption  is  actually  con¬ 
servative  except  in  the  downstream  transition 
zone,  and  hence  has  been  deemed  acceptable. 
Zufelt  (1999)  presents  a  test  that  can  be  used  to 
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determine  when  the  assumption  of  steady  state 
will  be  violated  to  a  sufficient  degree  that  unsteady 
methods  should  be  utilized.  Preliminary  results 
indicate  that  the  steady-state  assumption  is  accept¬ 
able  in  many  cases. 

A  number  of  expressions  are  available  for  cal¬ 
culating  the  equilibrium  ice  jam  thickness  in 
breakup  ice  jams,  that  is,  the  uniform  thickness 
section  of  an  ice  jam  that  has  formed  through 
shoving  of  the  ice  (Pariset  et  al.  1966;  Uzuner  and 
Kennedy  1974;  Beltaos  1978, 1995).  The  expression 
selected  for  use  in  HEC-RAS  is  based  on 
the  ice-jam  force  balance  equation  in  which  the 
stresses  within  the  jam  are  balanced  against  the 
external  stresses  on  the  jam  (USACE  1997,  Daly 
et  al.  1998): 


0  =  ft£Pi 


(gpiSf-2Ci)Ti~TB  (17) 


where  |i  is  the  coefficient  of  internal  strength  of 
the  ice  cover,  sometimes  called  the  friction  coeffi¬ 
cient,  and  Cj  is  the  cohesion.  Cohesion  is  often 
neglected  for  breakup  ice  jams,  but  can  be  impor¬ 
tant  for  freezeup  jams.  The  shear  stress  on  the 
underside  of  the  jam,  x,  can  be  expressed  as 


x  =  pgRiSf. 


(18) 


For  negligible  cohesion  (Q  =  0),  eq  17  yields  the 
following  expression  for  the  equilibrium  ice  jam 
thickness: 


BSf 

T1=2n(l-Si) 


1  + 


I  |  4^i(l-s,)R1 

i  siBsi 


(19) 


+  ^|I1  =  PiSSfil  +  'C  (16) 

where  ax  is  the  stress  in  the  streamwise  direction 
x  (assumed  to  be  constant  across  the  width  of  the 
channel,  B),  r\  is  the  thickness  of  the  ice  accumula¬ 
tion,  xb  is  the  shear  resistance  provided  by  the  river 
banks,  and  x  is  the  shear  stress  on  the  underside 
of  the  jam.  This  force  balance,  expressed  by  Ashton 
(1986)  in  the  form  below,  was  used  by  Tuthill  et  al. 
(1998)  in  the  ICETHK  option  of  HEC-2: 


where  s j  is  specific  gravity  of  ice  (pj/p).  The  water 
surface  elevation  in  the  equilibrium  section  of  a 
breakup  ice  jam,  h,  is  calculated  by 

=  3/i  +  siTl*  (2°) 

In  summary,  application  of  equilibrium  ice 
thickness  models  (e.g.,  eq  19)  will  require  some 
knowledge  or  estimation  of  the  following  ice  prop¬ 
erties:  Q,  ji,  p,  pj,  nb,  tiy  nc,  L,  and  H,  of  which  two, 
p  and  pj,  are  well  known. 
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ICE  PROPERTIES  USED  IN  HYDRAULIC 
MODELS  OF  ICE  JAMS 

The  properties  of  ice  jams  included  in  the  mod¬ 
els  cited  above  include  ice  cover  roughness,  cohe¬ 
sion,  measures  of  internal  strength  such  as  inter¬ 
nal  friction  angle  and  the  angles  of  static  and 
dynamic  friction,  the  shear  force  on  the  ice  cover 
and  its  shear  strength,  and  frazil  particle  diameter 
and  fall  velocity.  Each  of  these  is  addressed  in  a 
separate  section  below,  followed  by  discussions 
of  ice  transport,  thermal  effects  on  ice  jams,  and 
miscellaneous  hydraulic  properties  of  ice  jams 
such  as  porosity,  anchor  ice  growth,  jam  erosion, 
ice  cover  breakup  and  movement,  thermal  effects, 
coefficient  of  ice  loss,  and  permeability. 

Ice  cover  roughness  (jjj) 

Ice  cover  roughness  is  the  most  thoroughly  ex¬ 
plored  hydraulic  property  of  ice  covers.  In  fact, 
the  biennial  River  Ice  Workshops  sponsored  by  the 
Canadian  government  (see,  e.g.,  Ismail  1997),  be¬ 
gan  as  a  workshop  on  the  hydraulic  resistance  of 
river  ice  in  1980  largely  devoted  to  developing  an 
understanding  of  ice  cover  roughness  (Tsang  and 
Beltaos  1980).  Yet  there  is  still  a  need  to  obtain  re¬ 
liable  field  measurements  of  ice  jams  roughness. 
A  review  of  formulas  for  computing  composite 
roughness  coefficients  as  well  as  a  discussion  of 
roughness  estimates  made  by  others  (Lotter  1932, 
1941;  Belokon  1938, 1940;  Panov  1960)  is  included 
in  Nezhikhovskiy  (1964).  He  suggests  that  many 
published  values  of  n j  and  nc  are  incorrect  due  to 
errors  in  measured  and  estimated  discharge  and 
water  slope.  Also,  the  smoothing  of  the  bottom 


surface  of  the  ice  over  time  and  as  the  ice  cover 
thickens — effects  noted  by  many  observers — can 
be  a  source  of  error  in  estimates  of  roughness, 
which  is  often  taken  as  constant  over  time.  Al¬ 
though  it  has  been  argued  by  some  that  excessive 
reliance  is  placed  on  this  single  work  by 
Nezhikhovskiy,  a  careful  review  of  this  work  in¬ 
dicates  that  both  the  method  and  the  number  of 
samples  utilized  by  Nezhikhovskiy  inspire  confi¬ 
dence.  Rather,  it  is  the  incorrect  application  of  re¬ 
sults  clearly  intended  by  Nezhikhovskiy  for  use 
during  initial  ice  cover  formation  due  to  frazil  or 
sheet  ice  accumulation  or  frazil  deposition  only, 
to  other  situations,  particularly  breakup  jams,  that 
engenders  difficulty  in  roughness  estimation. 

Roughness  is  approached  in  several  ways  and 
using  several  coefficients.  Roughness  is  almost 
always  determined  using  a  two-layer  system,  that 
is,  dividing  the  flow  area  into  separate  layers  af¬ 
fected  by  the  bed  and  the  ice  cover  (Fig.  5).  A  com¬ 
posite  channel  roughness  coefficient  is  calculated 
through  application  of  standard  flow  equations 
such  as  Manning's  equation  or  the  Darcy- 
Weisbach  equation  (numerous  researchers)  or 
boundary-layer  theory  (e.g.,  Gogus  and  Tatinclaux 
1981).  Given  an  estimate  of  bed  roughness,  the 
roughness  of  the  underside  of  the  ice  can  then  be 
calculated  using  a  formula  such  as  the  Belokon- 
Sabaneev  formula  (eq  11).  Roughness  has  also  been 
calculated  from  estimates  of  the  shear  stress  on 
the  underside  of  the  ice  cover  (e.g.,  Knowles  and 
Hodgins  1980,  Hara  et  al.  1997).  Unfortunately,  few 
data  useful  in  calculating  roughness  exist,  particu¬ 
larly  for  breakup  jams.  Also,  the  presence  of  an 
ice  jam  can  cause  changes  in  bed  roughness,  which 
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in  turn  affects  the  calculation  of  ice  roughness.  In 
practice,  then,  roughness  is  generally  estimated 
during  calibration  of  numerical  models,  where  it 
essentially  becomes  a  lumped  variable  represent¬ 
ing  other  unknown  jam  components. 

The  most  commonly  used  roughness  coefficient 
is  the  Manning's  n,  followed  by  the  Darcy  friction 
factor  /.  The  Chezy  C  and  the  equivalent  sand 
roughness  k  are  occasionally  used  as  well,  but  less 
so  in  recent  years,  and  they  are  not  included  in 
the  most  frequently  used  hydraulic  models.  The 
following  discussion  will  be  limited  to  use  of  the 
various  Manning's  n  values  for  the  underside  of 
the  ice,  nx,  and  the  composite  roughness,  nc  Nu¬ 
merous  methods  of  determining  nc  are  given  in 
Uzuner  (1975).  No  attempt  has  been  made  to  con¬ 
vert  roughness  estimates  reported  in  the  form  of 
/,  C,  or  k  to  n.  Table  1  and  Figure  6  summarize  re¬ 
ported  values  of  nY  and  nc. 

Roughness  calculated  from  discharge 
measurements  or  estimates 

Based  on  about  500  observations  between  1936 
and  1959,  Nezhikhovskiy  reports  that  ni  for 
smooth  ice  covers  (i.e.,  generally  formed  through 
heat  transfer  processes  as  opposed  to  accumula¬ 
tion)  ranges  from  0.010  to  0.012  during  initial  ice 
cover  formation  and  from  0.008  to  0.010  thereaf¬ 
ter.  A  total  of  368  measurements  at  19  stations  on 


14  Russian  rivers  meeting  the  following  criteria 
was  used  to  determine  n{  for  rough  ice  covers:  the 
ice  cover  was  either  rough  or  had  frazil  deposits, 
sufficient  measurements  were  made  each  year  to 
determine  time-related  effects,  and  slope  was 
measured.  The  following  equation  obtained  from 
the  Manning  formula,  assuming  =  yx/2,  was 
used  to  calculate  nc : 


nc  =  0.63 


(21) 


and  ttj  was  then  back-calculated  using  eq  13. 

Based  on  his  results,  Nezhikhovskiy  concluded 
that  the  roughest  ice  covers  were  observed  dur¬ 
ing  the  first  three  days  of  ice  cover  and  that  n-x  tends 
to  decrease  over  the  course  of  the  winter.  Loose 
frazil  deposits  tended  to  smooth  more  rapidly  than 
dense  deposits.  He  also  noted  that  the  bottom  of 
frazil  deposits  tended  to  be  wavy,  and  that  thicker 
deposits  had  larger  amplitude  ridges;  both  thick 
and  thin  deposits  smoothed  during  the  winter.  By 
categorizing  ice  covers  into  three  categories  (ice 
covers  formed  primarily  from  loose  frazil,  dense 
frazil,  and  sheet  ice),  Nezhikhovskiy  developed 
the  curves  shown  in  Figure  7  to  describe  n-x  at  the 
time  of  initial  ice  cover  formation  in  terms  of  thick¬ 
ness.  In  the  case  of  frazil  accumulations  or  depos- 


Table  1.  Values  of  ice  roughness  coefficient  (iij)  and  composite  roughness  coefficient  ( nc ) 
calculated  from  discharge  measurements. 


(n  i) 

fnc) 

Comment 

Reference 

0.010-0.012 

Sheet  ice,  early  winter* 

Nezhikhovskiy  (1964) 

0.008-0.010 

Sheet  ice,  late  winter 

Nezhikhovskiy  (1964) 

0.010-0.06+ 

Ice  cover  formed  from  loose  frazil* 

Nezhikhovskiy  (1964) 

0.013-0.09+ 

Ice  cover  formed  from  dense  frazil* 

Nezhikhovskiy  (1964) 

0.015-0.10+ 

Ice  cover  formed  from  sheet  ice* 

Nezhikhovskiy  (1964) 

0.010-0.028** 

0.018-0.027 

Sheet  ice 

Carey  (1966) 

0.004-0.013** 

0.015-0.022 

Sheet  ice 

Carey  (1967) 

0.10 

0.090-0.109 

Breakup  jams 

Beltaos  (1978) 

0.057-0.065, 

0.041-0.046 

Breakup  jam 

Andres  (1980) 

Hj  =  0.060 

0.010-0.015 

Breakup  jam 

Knowles  and  Hodgins  (1980) 

0.053-0.142 

Breakup  jams 

Michel  (1980) 

0.013-0.040 

Freezeup  jam* 

Beltaos  (1981) 

0.033-0.041++ 

Freezeup  jam* 

Beltaos  (1983) 

0.072 

Breakup  jams 

Andres  and  Doyle  (1984) 

0.020-0.15 

Freezeup  jam,  frazil  deposits 

Majewski  and  Grzes  (1986) 

*Within  three  days  of  formation. 

+Higher  values  for  thicker  accumulations. 
**Lower  values  earlier  in  the  winter, 
ft  Higher  values  for  thinner  accumulations. 
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Korbaylo  and  Shumilak  (1999) 
White  and  Acone  (1 998) 
Tuthil!  and  White  (1997) 
Tuthill  et  a!.  (1996) 
McGilvary  et  al.  (1995) 
Rivard  et  a!.  (1984) 
Vogel  and  Root  (1981) 
Andres  and  Doyle  (1984) 
Knowles  and  Hodgins  (1980) 
Andres  (1980) 
Beltaos  (1978) 
Nezhikhovskiy  (1964) 
Tuthill  and  Mamone  (1998) 
Majewski  and  Grzes  (1986) 
Beltaos  (1983) 
Beltaos  (1981) 
Nezhikhovskiy  (1964) 
Nezhikhovskiy  (1964) 


0  0.02  0.04  0.06  0.08  0.10  0.12  0.14  0.16 


Ice  Roughness,  /7- 


Figure  6.  Reported  values  of  ice  roughness. 


0.1  1  10  (m) 

Ice  Thickness 


Figure  7.  Nezhikhovskiy 's  (1964)  suggested  coefficients  of  ice  roughness  (r\j)  at  the 
time  of  initial  ice  cover  formation,  based  on  the  thickness  and  type  of  the  parent  ice 
material 
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Determined 
By  Modeling 


its,  thickness  is  taken  as  the  measured  or  estimated 
thickness  of  the  ice  cover,  whereas  for  accumula¬ 
tions  made  up  of  broken  sheet  ice,  thickness  is  that 
of  the  parent  ice  material  He  specifically  states 
that  these  data  are  applicable  to  initial  ice  cover 
formation  and  not  to  breakup  jams.  Because 
the  data  underlying  this  figure  were  carefully 
collected  and  the  sample  size  was  large  (n  =  368), 
they  provide  a  valuable  resource  when  selecting 
roughness  values  at  the  time  of  initial  ice  cover 
formation. 

Beltaos  (1981)  calculated  nx  for  a  January  1980 
frazil  ice  jam  on  the  Thames  River,  New 
Brunswick,  using  estimated  discharge  and  thick¬ 
ness.  Like  Nezhikhovskiy,  he  observed  that  the 
initial  roughness  was  the  highest  and  that  smooth¬ 
ing  occurred  over  time,  also  three  days  in  this  case. 
Although  four  methods  of  estimating  thickness 
and  roughness  were  presented,  two  were  consid¬ 
ered  plausible.  The  values  of  nx  ranged  from  0.040 
to  0.016  using  equilibrium  ice  jam  criteria  and  from 
0.038  to  0.013  using  Nezhikhovskiy 's  relationship 
between  and  thickness.  At  the  time  of  initial 
freezeup  in  1981,  Beltaos  (1983)  estimated  that 
roughness  was  between  0.033  and  0.041,  depend¬ 
ing  on  thickness. 

Measurements  on  the  Vistula  River,  Poland, 
during  January  1982  allowed  Majewski  and  Grzes 
(1986)  to  determine  roughness  values  for  a  frazil 
ice  jam.  Over  the  course  of  20  days,  they  found 
0.020  <  rii  <  0.15,  with  the  roughest  sections  asso¬ 
ciated  with  shoved  areas.  Estimates  based  on  field 
measurements  made  in  1984  and  1985  ranged  from 
fti  =  0.012  for  sheet  ice  to  n *  =  0.080  for  frozen  ice 
floes  mixed  with  frazil.  They  also  noted  the  de¬ 
crease  in  roughness  over  time. 

Carey  (1966,  1967)  made  careful  observations 
of  the  sheet  ice  cover  near  a  USGS  gage  on  the  St. 
Croix  River,  Wisconsin,  during  three  winters 
(1964-66).  He  noted  dunes  and  ripples  on  the  un¬ 
derside  of  the  ice,  and  observed  that  the  wave¬ 
length  and  amplitude  of  the  dunes  and  ripples 
near  the  bank  were  greater  than  for  ice  taken  from 
midstream.  In  contrast  to  Nezhikhovskiy's  work 
with  frazil  deposits  and  accumulations  of  frazil 
and  sheet  ice,  Carey  found  that  roughness  for  sheet 
ice  covers  tends  to  increase  after  ice  cover  forma¬ 
tion.  Using  his  own  formulation  for  calculating  n-x 
and  nc,  he  reported  n j  between  0.010  and  0.028  and 
nc  between  0.018  and  0.027  during  1964  and  1965, 
with  lower  values  reported  earlier  in  the  winter 
(Carey  1966).  Later,  Carey  (1967)  indicated  that 
these  values  may  have  been  underreported  be¬ 
cause  of  an  incorrect  assumption  in  his  method. 


Using  a  modified  approach,  he  computed  0.004  < 
n-x  <  0.013  and  0.015  <  nc  <  0.022  for  measurements 
made  during  1966. 

Michel  (1980)  states  that  the  roughness  of 
breakup  jams  will  be  primarily  a  function  of  the 
shape  and  size  of  the  floes  making  up  the  jam.  For 
example,  he  points  out  that  jams  made  up  largely 
of  small  floes  will  have  smaller  roughness  than 
those  made  up  of  large  floes.  Although  Michel  dis¬ 
putes  the  use  of  the  relationship  suggested  by 
Nezhikhovskiy,  in  fact  his  reliance  on  floe  shape 
and  size  is  related  to  Nezhikhovskiy's  reliance  on 
parent  ice  thickness.  Michel  suggests  that  in  the 
case  of  breakup  jams  that  are  decaying,  the  rough¬ 
ness  will  decrease  with  decay.  Using  his  own  for¬ 
mulation,  he  determined  values  of  nc  for  ten  loca¬ 
tions  on  three  Canadian  Rivers  ranging  from  0.053 
to  0.142. 

Using  an  approach  combining  elements  of  the 
Nezhikhovskiy  approach,  Beltaos  (1978)  deter¬ 
mined  nc  and  n j  for  two  breakup  ice  jams  on  the 
Smoky  and  Wapiti  Rivers,  Alberta.  The  ice  rough¬ 
ness,  ttj,  was  found  to  be  constant  and  equal  to 
0.10,  while  nc  varied  from  0.090  to  0.109.  Using  a 
similar  approach,  Andres  (1980)  estimated  0.057 
<nx<  0.065  {ni  =  0.060)  and  0.041  <nc<  0.046  for 
the  April  1978  Athabasca  River  (Alberta)  breakup 
ice  jam  based  on  measured  discharge,  slope,  and 
stage.  When  data  from  this  jam  were  combined 
with  1977  and  1978  data,  Andres  and  Doyle  (1984) 
calculated  an  average  nx  =  0.072. 

The  shear  stress  approach  was  used  by  Knowles 
and  Hodgins  (1980)  to  estimate  n j  and  the  Belokon- 
Sabaneev  equation  to  calculate  nc  for  two  breakup 
jams  on  the  Thames  River,  New  Brunswick.  They 
found  that  nx  ranged  from  0.01  to  0.015  and  de¬ 
pended  on  the  bed  roughness  and  the  ratio  of  jam 
thickness  to  depth. 

Roughness  estimated  during  numerical 
model  calibration 

Vogel  and  Root  (1981)  modeled  breakup  jams 
on  the  Missisquoi  River  (Vermont)  using  nx  =  0.057. 
Rivard  et  al.  (1984)  used  HEC-2  to  model  a  May 
1983  breakup  ice  jam  on  the  Mackenzie  River.  The 
roughness  (n-x)  of  the  intact  ice  sheet  was  assumed 
to  be  0.020.  The  best  fit  with  measured  water  sur¬ 
face  profiles  was  found  with  nx  =  0.045  within  the 
jam.  Similarly,  roughness  values  between  0.044 
and  0.090  provided  a  good  fit  in  modeling  ob¬ 
served  jam  stages  on  the  Aroostook  River  (White 
and  Acone  1998).  Tuthill  and  White  (1997)  used  n-x 
=  0.060  to  model  a  breakup  jam  on  the  Salmon 
River,  Connecticut.  The  same  value  was  report- 
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edly  used  by  Korbaylo  and  Shumilak  (1999)  in 
modeling  the  1986  Churchill  River  breakup  jam. 
Tuthill  and  Mamone  (1998)  used  a  lower  value,  n-x 
=  0.030,  when  modeling  freezeup  jams  on  the 
middle  Mississippi  River. 

Uncertainty  in  ice  roughness 

White  and  Daly  (1997)  examined  the  effects  of 
uncertainties  in  the  estimate  of  ice  roughness  on 


the  calculated  equilibrium  ice  jam  thickness.  They 
obtained  46  values  of  n±  from  calibrated  numeri¬ 
cal  models  of  three  breakup  ice  jams.  These  ranged 
from  very  smooth  (0.020)  to  very  rough  (0.150), 
with  a  mean  of  0.066  and  standard  deviation  of 
0.023  (Figure  8a).  The  distribution  was  found  to 
be  nearly  normal.  A  Monte-Carlo  simulation  of 
sample  size  10,000,  based  on  the  mean  of  0.066  and 
standard  deviation  of  0.023,  is  shown  in  Figure  8b. 


Breakup  Jam  “n"  Values 

a.  Distribution  of  Manning's  n  values  for  the  underside  of  the  ice  from 
calibrated  numerical  models  of  breakup  ice  jams  at  three  sites. 


Simulated  Breakup  Jam  “n"  Values 

b.  Distribution  of  simulated  Manning's  n  values.  (From  White  and  Daly  1997.) 
Figure  8.  Distribution  of  actual  and  simulated  Manning's  n  values. 
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The  95%  confidence  interval  is  (0.029, 0.105),  which 
can  be  interpreted  to  mean  that  one  can  be  95% 
confident  that  the  Manning's  n  values  for  breakup 
jams  will  lie  between  0.029  and  0.105,  based  on 
the  data  shown  in  Figure  8a. 

Using  these  data,  they  simulated  equilibrium 
ice  thickness  using  eq  19  for  a  rectangular  chan¬ 
nel  width  of  100  m  with  nb  of  0.030,  discharge  100 
m3s-1,  water  slope  0.001  m/m,  and  equilibrium 
jam  section  porosity  of  50%.  The  results  showed 
that  ice  thickness  is  relatively  insensitive  to 
changes  in  roughness  (i.e.,  large  differences  in 
roughness  resulted  in  relatively  small  differences 
in  ice  thickness).  Stage  was  found  to  be  more  than 
twice  as  sensitive  to  roughness  than  ice  thickness 
for  the  channel  geometry  tested.  However,  less  un¬ 
certainty  was  associated  with  stage  or  thickness 
than  with  roughness.  The  opposite  was  reported 
by  Healy  et  al.  (1997),  who  evaluated  the  non-equi¬ 
librium  breakup  jam  that  occurred  in  1986  on 
the  Thames  River  (New  Brunswick)  using  the  Ca¬ 
nadian  model  ICEJAM  (Fla to  and  Gerard  1986). 
They  found  that  thickness  was  more  sensitive  to 
changes  in  roughness  than  stage  in  this  situation: 
a  change  in  composite  roughness  from  0.040  to 
0.080  caused  65%  change  in  thickness  and  25% 
change  in  stage. 

Cohesion  (Q) 

Fish  and  Zaretsky  (1997)  state  that  "Ice  cohe¬ 
sion  is  the  principal  strength  parameter  of  ice." 
This  is  immediately  obvious  to  the  observer  of 
freezeup  jam  formation  and  progression,  in  which 
the  apparent  cohesive  strength  of  the  frazil  floes 
plays  a  large  part  in  ice  cover  accumulation.  In 
addition,  because  freezeup  jams  most  often  occur 
at  very  low  air  temperatures,  freezing  at  the  wa¬ 
ter  surface  level  within  the  accumulation  can  lead 
to  rapid  and  large  increases  in  strength.  Tuthill  et 
al.  (1998)  suggest  that  cohesion  can  range  up  to 
about  960  Pa  to  1200  Pa  in  freezeup  jams. 

Yet  cohesion  is  often  neglected,  particularly  for 
breakup  jams.  This  is  partly  because  no  actual  mea¬ 
surements  have  been  made.  In  both  ICETHK  and 
HEC-RAS,  Q  is  assumed  to  be  0  for  breakup  ice 
jams  formed  from  unconsolidated  rubble.  How¬ 
ever,  some  researchers  have  considered  cohesion 
in  breakup  jams.  Beltaos  (1978)  used  a  value  of  96 
Pa  for  cohesion  when  modeling  the  1976  breakup 
jams  on  the  Smoky  and  Wapiti  Rivers  (Alberta). 
Following  his  method,  Andres  and  Doyle  (1984) 
also  accounted  for  cohesion  in  their  analysis  of 
breakup  jams  on  the  Athabasca  River  (Alberta)  in 
1977, 1978,  and  1979,  assuming  the  same  value  of 


Q  =  96  Pa.  Assuming  a  value  of  p  =  1.3,  Calkins 
(1983)  calculated  Q  =  168  Pa  and  254  Pa  for 
breakup  jams  on  the  Ottauquechee  River  and  First 
Branch  White  River  (Vermont),  respectively.  He 
concluded  that  cohesion  in  breakup  jams  is  insig¬ 
nificant  compared  to  water  slope,  river  width,  and 
bankfull  depth.  Ashton  (1986)  reported  that  labo¬ 
ratory  results  suggest  that  Q  =  100  Pa  to  500  Pa. 

Coefficient  of  internal  strength  or 
friction  (p)  and  angle  of  internal  friction  (<|>) 

The  coefficient  of  internal  strength  of  the  ice 
cover,  or  p,  is  sometimes  called  the  friction  coeffi¬ 
cient.  This  coefficient  is  a  lumped  variable  repre¬ 
senting  the  effects  of  a  number  of  unknown  or  dif- 
ficult-to-determine  factors,  including  lateral  stress 
within  the  ice  cover,  ice  on  bank  friction,  angle  of 
internal  resistance,  porosity,  and  possibly  cohe¬ 
sion.  Michel  (1984)  suggests  that  p  =  1.3  and  <])  = 
30  should  be  used  when  modeling  breakup  jams, 
but  a  variety  of  values  have  been  determined 
(Table  2).  Based  on  observations  of  a  number  of 
jams,  Beltaos  (1981)  hypothesized  p  is  important 
in  modeling  thicker  jams,  but  less  important  for 
thinner  jams.  In  an  evaluation  of  11  case  studies 
of  breakup  jams,  Beltaos  (1982)  reports  a  range  of 
0.8  <  p  <  1.3,  with  an  average  of  1.2.  Two  other 
jams,  for  which  p  =  0.6  and  p  =  3.5,  were  presented 
but  not  included;  the  lower  value  was  based  on 
highly  uncertain  post-jam  observations  and  the 
upper  on  a  thin  jam  for  which  ice  on  bank  friction 
may  have  been  highly  influential.  Beltaos  (1983, 
1978)  calculated  p  =  1.6  for  a  1981  breakup  jam  on 
the  Thames  River,  New  Brunswick,  and  0.9  to  2.2 
for  breakup  jams  on  the  Smoky  and  Wapiti  Riv¬ 
ers,  Alberta.  He  used  a  value  of  p  =  1.2  in  model¬ 
ing  the  1992  and  1996  Credit  River  (Ontario) 
breakup  jams  (Beltaos  1999). 


Table  2.  Values  of  the  coefficient  of  internal 
strength  (p)  calculated  by  various  researchers. 


Range  of  \i 

Mean  value 

Reference 

0.9-2. 2 

Beltaos  (1978) 

1.0 

Andres  (1980) 

0.8-1.3 

1.2 

Beltaos  (1982) 

1.6 

Beltaos  (1983) 

1.6,  2.0 

Rivard  et  al.  (1984) 

1.06 

Prowse  (1986) 

0.8-2.7 

1.6 

Andres  and  Doyle  (1984) 

1.2 

Beltaos  et  al.  (1996) 

1.17 

Tuthill  and  White  (1997) 

1.5 

Korbaylo  and  Shumilak  (1999) 
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Andres  (1980)  calculated  p  =  1.0  for  the  1978 
breakup  jam  on  the  Athabasca  River  (Alberta). 
Later,  Andres  and  Doyle  (1984)  used  Beltaos' 
(1978)  approach  to  determine  breakup  jam  char¬ 
acteristics  for  the  Athabasca  River  breakup  jams 
in  1977  and  1979  as  well  as  1978,  reported  previ¬ 
ously.  Assuming  Q  to  be  96  Pa,  they  found  p  to 
range  from  0.8  to  2.7,  with  an  average  of  1.6.  How¬ 
ever,  they  point  out  that  p  is  very  sensitive  to  er¬ 
rors  in  discharge,  and  that  the  unexpectedly  wide 
range  of  p  is  probably  attributable  to  errors  in  es¬ 
timating  discharge.  Tuthill  and  White  (1997)  re¬ 
ported  p  =  1.17  when  modeling  a  breakup  jam  on 
the  Salmon  River,  Connecticut.  Prowse  (1986)  cal¬ 
culated  p  =  1.06  for  the  1983  ice  jam  at  the 
confluence  of  the  Liard  and  Mackenzie  Rivers, 
Northwest  Territories,  while  Rivard  et  al.  (1984) 
calculated  higher  values  (2.0,  1.6)  for  two  cross 
sections  of  the  Mackenzie  River  jam  that  same  year. 
Korbaylo  and  Shumilak  (1999)  used  a  value  of  p  = 
1.5  when  modeling  the  1986  Churchill  River 
breakup  jam.  Reported  values  of  p  are  summa¬ 
rized  in  Figure  9. 

If  the  lumped  variable  approach  is  not  used, 
the  friction  coefficient  is  calculated  using  the  fol¬ 
lowing  equation  using  theory  developed  for 
granular  materials: 

\i  =  k0Kp\(l-p)  (22) 


where  Atq  is  the  coefficient  of  friction  between  ice 
and  bank;  X  is  the  lateral  stress  coefficient,  some¬ 
times  called  /q;  and  Kp  is  the  passive  pressure  co¬ 
efficient  (ratio  of  longitudinal  to  vertical  stress  at 
the  time  of  failure),  labeled  kx  in  some  models.  The 
lateral  stress  coefficient  is  often  assumed  to  be  0.33 
(e.g.,  Tuthill  and  Mamone  1998).  The  coefficients 
k0  and  Kp  are  calculated  by 

*0  -  tan(<(>)  and  Kp  -  tan2  ^45  +  (23) 

where  <J)  is  the  angle  of  internal  friction  in  degrees. 
(J)  is  sometimes  approximated  using  the  angle  of 
repose  in  dry  materials,  and  often  assumed  to  be 
45°  in  ice.  Substituting,  eq  22  becomes 

(X  =  tan(<]>)  tan2  ^45  +  ^  jx,(l  -  p)  (24) 

in  which  the  dependence  of  the  friction  coefficient 
on  angle  of  internal  friction  can  be  seen. 

Michel  (1980)  calculated  an  average  value  of  <|> 
=  21°  based  on  data  from  10  jams.  Weiss  et  al.  (1981) 
report  11°  <  f  <  34°  in  a  series  of  shear  box  tests  on 
ice  rubble,  with  higher  values  for  rubble  formed 
from  thicker  parent  ice  sheets.  Slower  displace¬ 
ment  for  the  same  ice  thickness  also  resulted  in 
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Figure  9.  Calculated  values  of  the  coefficient  of  internal  strength  (p). 
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higher  <|).  Based  on  laboratory  shear  box  tests  on 
ice  rubble,  Prodanovic  (1979)  reports  <j)  =  53°  for 
thicker  ice  and  <j)  =  47°  for  thinner  ice.  Shear  box 
tests  by  Heilman  (1984),  cited  in  Ettema  and  Urroz- 
Aguirre  (1991),  reportedly  resulted  in  $  =  54°. 
Ettema  and  Urroz- Aguirre  (1991)  report  that  freeze 
bonding  can  affect  the  shear  and  compressive 
strengths  of  ice,  and  suggest  that  the  relatively  high 
values  of  ([)  often  reported  are  actually  an  effective 
friction  angle  reflecting  these  and  other  effects. 
Loset  and  Sayed  (1993)  conducted  shear  tests  in  a 
biaxial  compression  chamber  and  noted  that  sub¬ 
merged  freshwater  ice  rubble  froze,  resulting  in 
larger  pieces  and  higher  angles  of  internal  friction 
than  for  dry  rubble.  They  tested  different  size  dis¬ 
tributions  of  25-mm  and  100-mm  ice  blocks,  and 
found  that  uniform  size  distributions  produced 
higher  stresses  than  a  mixture  of  half  large  and 
half  small  blocks.  Timco  and  Cornett  (1999)  con¬ 
ducted  shear  testing  of  broken  saline  and  fresh¬ 
water  ice  rubble  using  a  biaxial  compression  cham¬ 
ber.  They  found  that  friction  angle  decreases  with 
increasing  strain  ratio  until  some  limiting  friction 
angle  is  reached.  For  freshwater  ice,  this  ranged 
from  about  76°  at  low  strain  ratio  (*1.2)  to  about 
38°  at  a  high  strain  ratio  (3.5).  They  hypothesized 
that  fracturing  of  the  rubble  at  the  higher  strain 
ratios  could  explain  the  change  in  friction  angle. 


Table  3.  Reported  values  of  the  angle  of  internal 
friction  (ty)  for  freshwater  ice. 


Range 
off l> 

Mean 

value 

Reference 

53° 

Prodanovic  (1979)  (thicker  ice) 

47° 

Prodanovic  (1979)  (thinner  ice) 

54° 

Heilman  (1984)* 

21° 

Michel  (1980) 

110-340 

Weiss  et  al.  (1981) 

0 

CO 

“p 

0 

3 

Loset  and  Sayed  (1993) 

38°-76° 

Timco  and  Cornett  (1999) 

*Cited  in  Ettema  and  Urroz- Aguirre  (1991). 


Reported  values  of  0  are  summarized  in  Table  3 
and  Figure  10. 

The  angle  of  internal  friction  and  the  passive 
pressure  coefficient  are  related  as  shown  in  eq  23 
(Fig.  11).  Some  researchers  report  Kp  rather  than 
()).  Beltaos  (1988)  has  calculated  values  of  Kp  =  8 
and  10  for  a  refrozen  breakup  jam  on  the  Thames 
River,  New  Brunswick,  with  Kp  =  9.6  providing 
the  best  model  results  using  RIVJAM.  Beltaos  and 
Burrell  (1991),  Beltaos  (1993),  and  Beltaos  et  al. 
(1996)  used  Kp  -  10,  12,  12,  and  12  when  using 
RIVJAM  to  model  four  other  breakup  jams.  For 
the  1992  and  1996  Credit  River  (Ontario)  breakup 
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Figure  10.  Reported  values  of  the  angle  of  internal  friction,  for  freshwater  ice. 
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Angle  of  Internal  Friction,  <f>  (*) 

Figure  11.  Relationship  between  angle  of  internal  friction,  and  passive  pressure 
coefficient,  Kp. 


jams,  Beltaos  (1999)  used  Kp  =  12.  Healy  et  al.  (1997) 
evaluated  the  non-equilibrium  breakup  jam  that 
occurred  in  1986  on  the  Thames  River  (New 
Brunswick)  using  the  Canadian  model  ICEJAM 
(Flato  and  Gerard  1986).  They  found  that  jam  thick¬ 
ness  and  stage  were  relatively  insensitive  to 
changes  in  Kp  as  it  was  varied  from  8  to  12.  Tuthill 
and  Mamone  (1998)  used  a  value  of  Kp  =  5.83,  cor¬ 
responding  to  ({)  —  45°,  when  modeling  freezeup 
jams  on  the  middle  Mississippi  River. 


variables.  In  some  cases,  density  (D)  of  frazil  de¬ 
posits  or  accumulations  has  been  reported.  Den¬ 
sity  and  porosity  are  related  as  follows: 


Reported  values  of  porosity  discussed  below  are 
summarized  in  Table  4  and  Figure  12. 


Porosity  (p) 

The  porosity  of  an  ice  accumulation  is  used  in 
estimating  the  volume  of  ice  contributing  to  the 
ice  accumulation  and  in  determining  such  vari¬ 
ables  as  its  coefficient  of  internal  strength.  Poros¬ 
ity  is  described  by 


where  Vv  is  the  volume  of  voids,  both  above  and 
beneath  the  water  surface,  V\  is  the  volume  of  the 
ice,  and  Vt  is  the  total  volume  of  the  jam.  As  Beltaos 
(1995)  noted,  the  porosity  of  an  ice  jam  can  be  much 
larger  than  the  porosity  of  the  floes  comprising 
the  jam.  Few  field  measurements  of  porosity  ex¬ 
ist;  most  reported  values  are  either  assumed  or 
back-calculated  from  measurements  of  other  ice 


Freezeup  ice  jams  and  frazil  deposits 

Petryk  (1995)  reports  that  ICESIM  uses  a  value 
of  0.73  in  modeling  porosity  of  ice  jams  formed 
from  undercover  deposition  of  frazil.  This  value 
was  in  the  range  of  0.65  to  0.85  that  Beltaos  (1995) 
calculated  from  data  reported  in  Kivisild  (1959). 
Beltaos  and  Dean  (1981)  made  bulk  measurements 
in  a  thick  (16-m)  frazil  deposit  and  found  porosi¬ 
ties  of  0.541  and  0.33  for  the  bottom  and  top  of  the 
deposit,  respectively.  They  also  reported  densities 
of  400  to  600  kg/m3,  with  lower  densities  near  the 
bottom  of  the  deposit  and  higher  densities  just 
below  the  surface  of  the  ice  cover.  Majewski  and 
Grzes  (1986)  measured  density  along  a  42-km-long 
freezeup  jam  on  the  Vistula  River  (Poland)  be¬ 
tween  WJbc/awek  and  P/ock  and  found  a  wide 
range  of  density  in  deposits  and  a  narrow  range 
of  density  in  shoved  areas.  They  reported  that 
density  of  hanging  dams  ranged  from  100  to  700 
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Table  4.  Reported  measured  or  estimated  values  of  porosity. 


Jam  type 

Porosity 

Reference 

Frazil  deposit 

0.33, 0.51 

Beltaos  and  Dean  (1981) 

(bulk  measurement) 

Ice  rubble 

0.19-0.50, 
p  =  0.34 

Weiss  et  al.  (1981) 

Frazil  deposit 

0.32-0.67 

Dean  (1986) 

Frazil  deposit 

0.24-0.89 

Majewski  and  Grzes  (1986) 

Frazil  accumulated 

0.35-0.45 

Majewski  and  Grzes  (1986) 

by  shoving 

Breakup  jam 

0.38  ±  0.1 

Prowse  (1990) 

Frazil  deposit 

p  =  0.59 

White  and  Lawson  (1992) 

(borehole  dilution  test) 

Frazil  deposit 

0.35-0.77, 

White  and  Lawson  (1992) 

(bulk  measurement) 

p  =  0.52 

Freezeup  jam 

0.40 

Shen  and  Wang  (1992) 

Frazil  deposits 

0.67  +  0.13 

Andersson  and  Daly  (1992) 

(laboratory) 

Breakup  jam 

0.70  ±  0.2* 

Lever  and  Gooch  (1998) 

(groundwater  flow  theory) 

*May  be  high  due  to  incomplete  grounding  or  piping  through  jam. 


kg/m3  (0.89  <  p  <  0.24),  while  the  shoved  sections 
of  the  jams  were  more  dense,  between  500  and  600 
kg/m3  (0.45  <p  <0.35). 

Although  Michel  (1984)  states  that  Kivisild's  re¬ 
ported  porosity  values  and  his  own  value  of  0.75 
represented  reasonable  values  for  freezeup  jams. 


smaller  values  of  frazil  deposit  porosity  are  re¬ 
ported  by  several  researchers.  Based  on  data  from 
borehole  dilution  tests  of  frazil  accumulations  in 
the  laboratory.  White  (1991)  estimated  mean  po¬ 
rosity  to  be  0.44,  using  the  Kozeny-Carmen  equa¬ 
tion  and  an  assumed  particle  diameter  of  0.2  mm. 
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Figure  12.  Reported  porosity  measurements  and  estimates . 
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Using  the  same  method.  White  and  Lawson  (1992) 
calculated  an  average  porosity  of  0.59  for  a  late- 
winter  frazil  deposit  on  the  Tanana  River,  Alaska. 
This  compared  favorably  to  porosities  determined 
by  direct  measurement  of  ice  and  water  volumes 
in  bulk  samples,  which  ranged  from  0.35  to  0.77, 
with  an  average  of  0.52.  Based  on  field  measure¬ 
ments,  Shen  and  Wang  (1992)  reported  porosity 
of  0.4  for  a  freezeup  jam  on  the  Yellow  River,  China 
(no  method  was  specified).  Tuthill  et  al.  (1998)  sug¬ 
gest  a  range  of  0.5  to  0.7  for  frazil  pans  at  freezeup 
when  using  ICETHK.  Tuthill  and  Mamone  (1998) 
used  p  =  0.5  when  modeling  a  freezeup  jam  on  the 
middle  Mississippi  River.  Measurements  of  the  po¬ 
rosity  of  frazil  ice  deposits  formed  on  trash  racks 
in  the  laboratory  were  made  by  Andersson  and 
Daly  (1992)  using  a  modified  CRREL  snow  den¬ 
sity  kit.  They  measured  20  samples  and  found  a 
mean  porosity  of  0.67  (Fig.  13). 

Breakup  ice  jams  and  ice  rubble 

Weiss  et  al.  (1981)  measured  porosity  for  manu¬ 
factured  ice  rubble  of  different  parent  ice  sheet 
thickness.  They  found  porosity  ranging  from  0.19 
to  0.50,  with  an  average  value  of  0.34.  Higher  po¬ 
rosity  was  generally  associated  with  thicker  par¬ 
ent  ice.  Prowse  (1988)  developed  a  method  to  esti¬ 


mate  porosity  of  breakup  ice  jams  based  on  the 
energy  and  mass  balance  of  an  ablating  ice  jam  at 
the  confluence  of  the  Liard  and  MacKenzie  Rivers 
(Northwest  Territories)  in  1983.  He  reported  an 
average  porosity  of  0.35.  Because  the  uncertainty 
of  the  method  is  largely  attributable  to  air  tem¬ 
perature,  he  varied  his  mean  temperature  by  ±  5°C, 
with  a  resulting  range  of  porosity  from  0.26 
(warmer  water  temperatures)  to  0.44  (cooler  wa¬ 
ter  temperatures).  Upon  improving  the  method 
(Prowse  and  Marsh  1989,  Prowse  1990),  the  range 
of  porosity  was  found  to  be  0.38  ±  10%,  quite  close 
to  the  commonly  assumed  value  of  0.4  (e.g., 
Beltaos  1978, 1995, 1999).  Michel  (1980, 1984)  and 
Tuthill  et  al.  (1998)  suggested  porosity  values  of 
0.4  to  0.5  when  modeling  breakup  jams.  Tuthill  et 
al.  (1996)  used  p  =  0.5  in  the  hydraulic  modeling 
of  a  breakup  ice  jam. 

Grounded  jams 

Porosity  is  particularly  important  in  the  case  of 
breakup  ice  jams  that  are  grounded,  because  per¬ 
meability  varies  with  porosity  and  hence  discharge 
will  also  vary  with  porosity,  affecting  backwater 
levels  upstream  from  a  grounded  jam.  Wong  et  al. 
(1985)  and  Beltaos  (1993, 1999)  used  a  theoretical 
model  of  porous  flow  presented  by  Bear  (1972)  in 
which  porosity  is  related  to  seepage  flow: 


Figure  13.  Porosity  measurements  of  frazil  deposits  formed  on  trashracks  in  the  CRREL  refrig - 
erated  flume.  (From  Andersson  and  Daly  1992.) 
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where  X  is  the  seepage  coefficient,  K  is  a  dimen¬ 
sionless  seepage  coefficient,  and  ds  is  a  variable 
describing  the  particles  ( ds  =  6/Ms,  where  Ms  is 
the  specific  area  per  unit  volume  of  the  solid  par¬ 
ticles,  or  the  ratio  of  surface  area  to  volume).  In 
tests  with  square  polyethylene  blocks,  Wong  et  al. 
(1985)  found  that  K  varied  from  0.3  to  0.7,  with 
the  higher  value  for  randomly  placed  blocks. 
Beltaos  (1999)  suggested  that  ds  can  be  approxi¬ 
mated  as  twice  the  average  ice  thickness  of  the  ice 
blocks.  This  approach  was  used  by  Lever  and 
Gooch  (1998),  who  obtained  careful  water  level 
measurements  and  discharge  estimates  at  rated 
cross  sections  for  a  breakup  jam  thought  to  be 
grounded  at  the  toe.  Assuming  K  =  0.6  and  ds  ~ 
twice  the  average  parent  ice  thickness,  they  calcu¬ 
lated  p  -  0.70  ±  0.2.  However,  this  high  value  for 
breakup  jam  porosity  may  be  a  result  of  incom¬ 
plete  grounding  or  piping  through  the  jam,  be¬ 
cause  the  analysis  includes  only  porous  flow. 

Shear  stress  on  the  underside  of  the  ice  cover  (x) 

As  noted  in  eq  18,  the  shear  stress  on  the  un¬ 
derside  of  the  ice  cover,  x,  is  a  function  of  p,  g,  Sf, 
and  the  variables  necessary  to  calculate  R{.  q,  nv 
nb,  and  nc.  Based  on  observed  water  profiles, 
Andres  (1980)  calculated  x  =  18.6  Pa  for  the  1978 
Athabasca  River  (Alberta)  breakup  jam.  Andres 
and  Doyle  (1984)  report  a  range  of  9.6  to  28.5  Pa 
and  an  average  of  17.9  Pa  for  shear  forces  experi¬ 
enced  during  breakup  jams  on  the  Athabasca  River 
(Alberta)  in  1977,  1978,  and  1979.  Rivard  et  al. 
(1984)  calculated  x  =  7.21  Pa  and  x  =  10.0  Pa  for 
two  cross  sections  of  the  1983  Mackenzie  River 
(Northwest  Territories)  ice  jam.  Dean  (1986)  esti¬ 
mated  shear  beneath  frazil  ice  deposits  to  be  about 
8  to  10  Pa  during  steady  midwinter  flow  and  60  to 
700  Pa  during  a  severe  water  level  change. 

Frazil  particle  diameter  (d)  and  fall 
velocity  coefficient  (Cf) 

Few  values  of  frazil  particle  sizes  exist  outside 
Alaska.  Gosink  and  Osterkamp  (1983)  reported 
frazil  particle  diameters  ranging  between  1  and  6 
mm  during  a  series  of  tests  on  the  Chatanika  River, 
Alaska.  Chacho  et  al.  (1986)  reported  frazil  diam¬ 
eters  of  5  mm  to  150  mm  in  the  Tanana  River.  White 
and  Lawson  (1992)  cited  unpublished  data  by 
Lawson  and  Brockett  that  report  particles  between 
1  mm  and  25  mm  in  diameter,  also  on  the  Tanana 
River,  Alaska.  They  found  mean  particle  sizes  be¬ 


tween  2  mm  and  8  mm  determined  by  sieve  analy¬ 
ses  on  the  samples.  Beltaos  and  Dean  (1981)  re¬ 
ported  smaller  diameter  particles  (1  mm  to  6  mm) 
in  a  frazil  deposit  on  the  Smoky  River,  Alberta, 
with  60%  by  weight  in  the  range  of  1.1  to  2.4  mm. 
Shen  and  Wang  (1992)  reported  mean  particle  di¬ 
ameter  of  about  1  cm  on  the  Yellow  River,  China. 
The  fall  velocity  coefficient  (Cf)  for  particles  5  mm 
to  1  cm  in  diameter  has  been  reported  to  be  about 
1.0  (Shen  and  Wang  1995). 


MISCELLANEOUS  VARIABLES  AFFECTING 
RIVER  ICE  HYDRAULICS 

Anchor  ice  growth 

Anchor  ice  growth  occurs  when  frazil  ice  ad¬ 
heres  to  rocks  or  other  materials  on  the  bed  of  a 
river.  The  presence  of  anchor  ice,  or  frazil  ice  ac¬ 
cumulated  on  submerged  surfaces,  can  affect  river 
hydraulics  in  a  number  of  ways,  through  increas¬ 
ing  the  effective  bed  level  and  changing  the  effec¬ 
tive  bed  roughness,  among  others.  Brown  et  al. 
(1953)  reported  that  frazil  ice  adhered  first  on  the 
upstream  face  of  stones  on  the  riverbed,  or  at  the 
point  of  attachment  to  the  bed  for  other  submerged 
objects.  Measurements  showed  that  the  average 
velocity  at  the  attachment  points  on  the  upstream 
face  of  the  stones  was  0.40  m/s  compared  to  0.77 
m/ s  along  the  sides  (where  anchor  ice  accumu¬ 
lated  later).  Under  favorable  conditions  the  anchor 
ice  grew  both  horizontally  and  vertically  to  form 
a  mat  of  anchor  ice.  Anchor  ice  thickness  appeared 
to  be  a  function  of  water  depth,  temperature,  and 
velocity,  and  reached  0.6  m  in  one  location. 

Data  from  Hirayama  et  al.  (1997)  from  the 
Niuppu  River,  Japan,  indicate  an  apparent  aver¬ 
age  density  of  about  0.5  g/cm3  for  anchor  ice.  They 
also  reported  that  the  density  of  anchor  ice  appears 
to  increase  slightly  as  water  velocity  increases. 
Using  a  six-hour  freezing  index  (FI,  units  °C  hr) 
for  the  hours  0000  to  0600,  they  measured  an  av¬ 
erage  anchor  ice  accumulation  rate  of  0.10  m3  per 
°C  hr.  Anchor  ice  accumulation  was  observed  to 
be  related  to  Froude  number,  F,  as  well  as  FI.  Ac¬ 
cumulations  occurred  for  0.2  <  F  <  1.5,  but  none 
were  seen  for  F  <  0.2.  A  high  FI  (i.e.,  more  intense 
cold)  allowed  frazil  ice  to  accumulate  at  lower 
values  of  F,  but  lower  FI  required  a  higher  F  be¬ 
fore  anchor  ice  accumulation  was  seen.  Finally, 
they  noted  that  more  anchor  ice  volume  results  on 
surfaces  with  higher  coefficients  of  static  friction. 

Kerr  et  al.  (1997,  1998)  investigated  anchor  ice 
growth  on  rounded  gravel  averaging  4.45  and 
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10.16  cm  in  diameter  in  a  refrigerated  flume.  They 
observed  a  definite  increase  in  effective  bed  eleva¬ 
tion  because  of  anchor  ice  formation,  but  the  pre¬ 
liminary  results  were  inconclusive  regarding 
roughness  effects.  Anchor  ice  was  observed  to  ini¬ 
tiate  and  accumulate  at  0.14  <  F  <  0.57.  Three  dif¬ 
ferent  types  of  anchor  ice  were  observed:  compact, 
relatively  stable  growth  in  the  form  of  "balls"  of 
frazil  attached  to  rocks  at  F  >  0.56;  a  loose,  rela¬ 
tively  unstable  growth  at  F  <  0.2  when  loose,  po¬ 
rous  tails  formed;  and  a  transitional  stage  during 
which  the  tails  either  flattened  out  to  form  a  mat 
or  were  released  (Fig.  14).  According  to  Daly,*  the 
loose  accumulations  were  very  unstable  and 
tended  to  release  in  large  sheets.  The  loose  accu¬ 
mulations  also  seemed  to  increase  effective  bed 
roughness  compared  to  the  compact  and  transi¬ 
tional-type  growth,  which  tended  to  smooth  the 
bed,  thereby  increasing  velocities  and  decreasing 
stage.  White  et  al.  (1999)  performed  a  series  of  tests 
in  the  same  flume  using  cobbles  of  median  diam¬ 
eter  4  cm.  They  reported  little  growth  at  a  Froude 
number  of  about  0.21  and  compact  or  transition- 
type  growth  that  eventually  formed  a  mat  of  an- 


*Personal  communication,  Steven  F.  Daly  Ice  Engineer¬ 
ing  Research  Division,  CRREL,  Hanover,  New  Hamp¬ 
shire,  1999. 


chor  ice  at  a  Froude  number  of  about  0.37,  within 
the  range  reported  by  Hirayama  et  al.  (1997). 

Jam  erosion 

Frazil  ice  and  ice  pieces  deposited  beneath  an 
ice  cover  can  erode  under  certain  conditions.  Ero¬ 
sion,  like  early  efforts  to  quantify  deposition  be¬ 
neath  ice  covers,  is  most  often  described  in  terms 
of  an  erosion  velocity.  Unlike  HEC-RAS,  ICETHK 
and  ICESIM  do  provide  for  erosion  of  deposited 
ice  through  the  use  of  a  critical  velocity  criterion. 
ICETHK  uses  the  following  equation,  developed 
from  continuity: 


H-\- 


(H-sp!) 


(28) 


where  T|t  is  the  ice  thickness  after  thinning  by  ero¬ 
sion  and  ve  is  the  critical  erosion  velocity  (Tuthill 
et  al.  1998).  The  critical  erosion  velocity  is  thought 
to  increase  over  the  course  of  a  winter  and  also  to 
vary  by  ice  type.  Tuthill  et  al.  (1998)  suggest  val¬ 
ues  of  vc  between  0.9  and  1.5  m/s  for  freezeup  jams 
and  1.2  <  ve  <  2.4  m/s  for  breakup  jams.  Petryk 
(1995)  reports  use  of  a  critical  erosion  velocity  of 
1.1  m/s  for  modeling  erosion  of  frazil  deposits  on 
the  Peace  River  using  the  S1MGLACE  model. 


Figure  14.  Anchor  ice  growth  characteristics  at  different  flow  velocities  and  flow  depths 
from  Kerr  et  al.  1998  (open  markers)  and  White  et  al  1999  (solid  markers). 
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Ice  cover  breakup  and  movement 

Well-documented  observations  of  ice  cover 
breakup,  movement,  and  jamming  are  provided 
by  Beltaos  (1981)  for  the  1980  Thames  River  (New 
Brunswick)  breakup.  The  candled,  relatively  weak 
ice  cover  was  observed  to  lift,  crack,  and  move 
downstream  at  a  rate  of  about  0.6  to  0.9  m/s. 
Downstream  movement  averaged  about  0.8  m/ s, 
with  some  consolidation  noted  at  0.3  m/s.  Once 
jammed,  the  jam  progressed  upstream  at  about  0.9 
m/s.  The  velocity  of  two  different  jam  releases  (in 
different  places)  was  also  about  0.9  m/s.  In  1981, 
the  ice  cover  was  observed  to  lift  and  crack,  but 
not  break,  for  water  level  rises  of  1  to  2  cm/hr,  3  to 
5  cm/hr,  5.5  cm/hr,  7  to  9  cm/hr,  and  8  cm/hr  at 
different  times  (Beltaos  1983).  Surface  movement 
of  broken  ice  was  observed  at  velocities  of  up  to 

1.2  m/s.  The  velocity  of  an  ice  run  caused  by  the 
failure  of  an  upstream  jam  was  about  1.7  m/s.  The 
same  velocity  was  reported  by  Hicks  et  al.  (1997) 
for  an  ice  run  caused  by  the  failure  of  a  1993  ice 
jam  on  the  St.  John  River,  Maine  and  New 
Brunswick. 

Andres  and  Doyle  (1984)  reported  average  sus¬ 
tained  velocity  of  the  ice  cover  breakup  front  on 
the  Athabasca  River  (Alberta)  to  be  3.7  m/s  and 

3.3  m/s  during  1978  and  1979,  respectively.  Higher 
values  were  observed  for  short  durations:  5.5  m/ 
s  in  1977  and  1979,  and  5.0  m/s  in  1978.  Observa¬ 
tions  of  the  speed  of  the  peak  discharge  and  the 
ice  front  following  the  failure  of  a  1995  breakup 
jam  on  the  Porcupine  River,  Yukon,  were  reported 
by  Jasek  (1997).  He  noted  that  the  peak  discharge 
arrived  faster  (5.4  m/s)  than  the  ice  front  (4.2  m/ 
s).  Ferrick  et  al.  (1996)  used  a  Doppler  radar  sys¬ 
tem  to  measure  the  movement  of  a  breakup  ice 
run  (March  1993)  and  a  frazil  ice  run  (December 
1993)  on  the  Connecticut  River,  New  Hampshire. 
Velocities  during  the  breakup  run  varied  from  1 
m/s  to  2.5  m/s  with  a  mean  velocity  of  1.35  m/s. 
The  frazil  run  was  slower,  ranging  from  0.5  m/s 
to  0.65  m/s  with  a  mean  of  0.54  m/ s,  which  com¬ 
pared  well  with  estimates  made  from  video  record¬ 
ings  that  averaged  0.52  m/s. 

Coefficient  of  ice  loss  (q) 

The  breakup  and  movement  of  an  ice  cover  of¬ 
ten  results  in  the  stranding  of  ice  along  the  banks 
in  the  form  of  shear  walls  or  ice  that  was  carried 
overbank  and  deposited  in  the  floodplain.  It  is  also 
likely  that  warm  temperatures  during  breakup 
may  melt  ice  pieces  travelling  relatively  long  dis¬ 
tances.  Ice  loss  due  to  stranding  and  melting  can 
decrease  the  ice  supply  available  for  downstream 


jamming  and  thus  is  of  interest  when  estimating 
ice  jam  volume.  Ice  loss  during  breakup  and  ice 
cover  movement  can  be  described  using  an  ice  loss 
coefficient,  q: 

Vt-Vj 

r,  — _ L  /r^rw 


where  Vt  is  the  total  upstream  ice  supply  and  Vj  is 
the  volume  of  ice  contained  in  the  jam. 

The  ice  loss  coefficient  is  expected  to  be  low  for 
shallow,  steep  rivers  with  short  reaches  of  contrib¬ 
uting  ice  cover  and  higher  for  more  gently  slop¬ 
ing  rivers  with  longer  reaches  of  contributing  ice 
cover,  many  tributaries,  and/or  low  banks 
(USACE  1982).  Using  a  porosity  of  0.40,  Prowse 
(1986)  estimated  that  loss  coefficients  ranged  from 
0.90  to  0.93  for  ice  jam  volume  and  ice  supply  data 
collected  by  Calkins  (1978)  on  the  relatively  short, 
steep  White  River  (Vermont).  Based  on  data  from 
Beltaos  (1983)  and  assuming  a  porosity  of  0.40, 
Prowse  (1986)  estimated  a  loss  coefficient  of  0.81 
for  the  Liard  River  mainstem  (British  Columbia 
and  Northwest  Territories)  during  the  breakup  of 
1983,  but  he  also  reported  an  ice  loss  of  about  40% 
in  one  day  as  ice  melted,  became  stranded,  or  was 
broken  into  small  enough  pieces  to  be  transported 
beneath  the  jam. 

Shear  strength 

Beltaos  and  Dean  (1981)  measured  shear 
strength  within  a  thick  (up  to  16  m)  frazil  deposit 
on  the  Smoky  River,  Alberta,  using  an  experimen¬ 
tal  shear  vane.  They  found  that  shear  strength 
varied  from  year  to  year  and  also  that  it  tends  to 
increase  with  height  above  the  bottom  of  the  ac¬ 
cumulation.  Reported  values  ranged  up  to  about 
80  kPa.  According  to  Dean  (1986),  the  shear 
strength  of  loose,  freshly  deposited  frazil  is  on  the 
order  of  10  Pa,  while  an  undisturbed,  undegraded 
frazil  accumulation  will  have  a  shear  strength  of 
about  30  to  60  kPa.  Compacted  frazil  accumula¬ 
tions  are  said  to  attain  a  shear  strength  of  about  75 
kPa.  Thermal  degradation  will  reduce  shear 
strength  in  frazil  deposits  to  between  1  and  10  kPa. 

Solid  ice  covers  have  much  higher  levels  of 
shear  strength  than  frazil  deposits.  Based  on  ob¬ 
servations  of  the  speed  of  the  moving  ice  front  for 
breakup  ice  jams  on  the  Athabasca  River  (Alberta) 
during  1977,  1978,  and  1979,  Andres  and  Doyle 
(1984)  estimated  the  strength  of  the  ice  cover  at 
between  310  and  650  kPa  for  1-m  nominal  ice  thick¬ 
ness.  They  reported  that  these  values  were  reason- 
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able  compared  to  estimates  made  by  others  from 
force  measurements  on  piers  at  a  downstream 
bridge. 

Thermal  effects  on  ice  jams 

Once  an  ice  cover  or  ice  accumulation  has 
formed,  it  will  thicken  largely  because  of  heat 
transfer  processes,  which  Ashton  (1986)  identified 
as  shortwave  (solar)  radiation,  longwave  radia¬ 
tion,  evaporation  or  condensation,  convection,  and 
precipitation.  Ice  cover  decay  is  also  a  function  of 
heat  transfer  processes,  but  it  is  more  complicated 
because  melting  at  both  the  top  and  bottom  sur¬ 
face  of  the  ice,  and  melting  within  the  ice  cover 
matrix  itself  are  involved.  Calculating  a  detailed 
heat  budget  that  describes  the  decay  process  can 
be  quite  complex,  because  it  depends  on  latitude, 
time  of  year,  time  of  day,  hours  of  sunlight,  cloud 
cover,  wind  speed,  air  temperature,  water  tem¬ 
perature,  humidity,  ice  type,  and  other  variables. 
Solar  radiation  plays  an  important  role  in  the 
heat  budget,  particularly  in  the  spring  when  the 
hours  of  daylight  are  increasing,  because  melting 
within  the  ice  cover  is  largely  caused  by  in¬ 
creased  absorption  of  solar  radiation.  This  inter¬ 
nal  melting  process  reduces  the  internal  strength 
of  the  ice  cover  and  its  resistance  to  failure  at  low 
stress  levels. 

Heat  flux 

Huokuna  (1988)  reported  on  the  modeling  of 
ice  cover  growth  and  decay  on  a  37.6-km  reach  of 
the  Oulujoki  River,  Finland.  The  heat  flux  between 
the  water  and  air  includes  contributions  from 
short-  and  longwave  radiation,  convection,  and 
evaporation.  A  heat  flux  of  2  W/m2  was  assumed 
between  the  water  and  the  channel  bed.  The  heat 
flux  between  the  water  and  the  ice  cover  is  a  func¬ 
tion  of  flow  turbulence  and  water  temperature. 
The  model  assumed  the  thermal  conductivity  of 
the  (snow-covered)  ice  cover  to  be  1.8  W/mK, 
compared  to  2.24  W/mK  for  pure  ice.  Calculated 
ice  thickness  resulting  from  thermal  decay  at  two 
cross  sections  compared  favorably  to  observed  ice 
thickness.  Careful  measurements  of  the  decay  of 
the  pre-breakup  ice  cover  and  the  1983  jam  at  the 
confluence  of  the  Liard  and  Mackenzie  Rivers 
(Northwest  Territories)  provided  information  used 
in  determining  the  magnitude  of  the  major  heat 
fluxes  contributing  to  ice  decay  (Prowse  1988, 1990; 
Prowse  and  Marsh  1989).  Between  6  May  and  9 
May  1983,  approximately  3.24  x  107  m3  of  ice 
melted. 


Albedo 

The  absorption  of  solar  radiation  by  an  ice  cover 
is  a  function  of  the  reflectivity  of  its  surface.  Dif¬ 
ferent  surfaces  (e.g.,  water,  snow,  ice)  reflect  solar 
radiation  to  different  degrees.  The  ratio  of  the  so¬ 
lar  radiation  reflected  by  a  given  surface  to  the 
incoming  solar  radiation  is  termed  the  albedo  of 
the  surface.  A  high  albedo  indicates  high 
reflectivity  (and  low  absorption),  while  a  low  al¬ 
bedo  indicates  low  reflectivity  (and  a  high  degree 
of  absorption).  The  albedo  depends  on  both  the 
condition  of  the  surface  (i.e.,  rough,  smooth)  and 
the  sun  angle.  According  to  Williams  and  Gold 
(1963),  the  albedo  of  new  snow  is  about  0.8  to  0.9, 
compared  to  an  albedo  of  0.4  to  0.6  for  melting 
snow,  and  0.05  to  0.15  for  water.  Gray  and  Male 
(1981)  report  albedos  of  0.75  to  0.95  for  fresh  snow, 
0.3  to  0.4  for  sea  ice,  and  0.05  to  0.30  for  water. 
Kondratyev  (1954),  as  presented  in  Gray  and  Male 
(1981),  reports  an  inverse  relationship  between  sun 
angle  and  albedo.  Measured  albedo  increased  from 
0.86  to  0.95  for  compact,  dry,  clean  snow  as  the 
sun  angle  decreased  from  30.3  to  25.1  degrees. 
They  reported  smaller  increases  for  wetter  and 
more  porous  snow.  Measurements  of  albedo  for 
brash  ice  reported  in  Prowse  and  Marsh  (1989) 
ranged  from  0.04  to  0.15,  with  a  mean  value  of  0.08. 
Candled  thermally  grown  ice  had  an  albedo  of 
0.39,  while  granular  thermally  grown  ice  had  an 
albedo  of  0.55.  During  an  ice  dusting  operation 
on  the  Platte  River,  Nebraska,  the  albedos  of 
dusted  and  non-dusted  areas  were  measured  be¬ 
tween  1  March  and  7  March  1979  (USAED  Omaha 
1979).  The  average  albedo  of  the  non-dusted  ice 
decreased  during  the  measurement  period  from 
about  0.65  on  1  March  to  0.40  on  5  March,  prob¬ 
ably  because  of  the  formation  of  pools  of  water  on 
the  ice  surface  from  rainfall,  snowmelt,  and  ice 
melting.  The  albedo  of  the  dusted  ice  was  about 
0.2  to  0.25. 

Photosynthetically  active  radiation  (PAR) 

The  presence  of  an  ice  cover  is  thought  to  pre¬ 
vent  or  minimize  photosynthesis  through  the  ab¬ 
sorption  or  reflection  of  solar  radiation  in  the  form 
of  visible  light.  The  wavelengths  from  approxi¬ 
mately  380  to  760  nm  are  used  by  photosynthetic 
organisms  located  throughout  the  water  column 
(Brock  et  al.  1984,  Atlas  and  Bartha  1998),  with  the 
range  between  400  and  700  nm  known  as  photo¬ 
synthetically  active  radiation  (PAR).  Bolsenga  et 
al.  (1991, 1996)  reported  that  transmittance  of  PAR 
through  clear  freshwater  ice  ranged  from  less  than 
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2.5%  for  snow-covered  ice  up  to  85%  through  21- 
cm-  (8.3-in.-)  thick  clear  ice.  Transmittance  through 
clear  ice  on  both  clear  and  cloudy  days  was  slightly 
higher  for  clear  days  in  the  upper  2  m  (6.6  ft)  of 
depth  and  was  not  significantly  different  below 
that  point.  No  actual  values  of  incident  or  trans¬ 
mitted  PAR  were  reported.  Vincent  et  al.  (1998) 
reported  transmittance  of  -1%  to  -10%  in  four 
Antarctic  freshwater  lakes,  with  PAR  of  about  80 
jimols/m2  at  10-m  (32.8-ft)  depth  in  Lake  Vanda 
and  about  8  |imols/m2  at  5-m  (16.4-ft)  depth  in 
Lake  Bonney.  White  and  Melloh  (1999)  measured 
PAR  at  mid-depth  in  a  shallow  ice-covered  river 
in  Vermont  and  found  that  the  mean  daily  inte¬ 
grated  irradiance  was  3.3  mols/m2/day  during  a 
late-winter  ice-covered  period.  Their  maximum 
observed  PAR  of  488  fimols/m2  is  comparable  to 
open-water  wintertime  values  reported  by  Guasch 
et  al.  (1998)  for  rivers  in  Spain  as  well  as  to  sum¬ 
mertime  measurements  for  shaded  rivers. 

Permeability 

Dean  (1976)  and  Beltaos  and  Dean  (1981)  used 
a  constant  head  permeameter  to  determine  the  in¬ 
trinsic  permeability  of  frazil  deposits.  Dean  re¬ 
ported  a  value  of  1.53  x  10“5  cm2,  and  Beltaos  and 
Dean  reported  1.56  x  10"5  ±  9.19  x  10~7  cm2.  These 
are  comparable  to  the  intrinsic  permeability  of  a 
coarse  sand  or  fine  gravel  deposit.  However,  the 
test  method  involved  destructive  sampling  that 
may  have  affected  their  results.  White  (1991)  used 


a  borehole  dilution  method  to  determine  perme¬ 
ability  of  a  frazil  ice  accumulation  in  the  labora¬ 
tory.  She  reports  an  intrinsic  permeability  of  6.1  ± 
1.94  x  10-7  cm2  (coefficient  of  variation  32%)  and 
hydraulic  conductivity  of  3.32  ±  1.07  x  10“2  cm/s 
(coefficient  of  variation  32%).  These  values  corre¬ 
spond  to  a  silty-to-clean  sand  deposit.  Similar  tests 
on  frazil  ice  deposits  in  the  Tanana  River,  Alaska, 
(White  and  Lawson  1992)  resulted  in  the  calcula¬ 
tion  of  an  intrinsic  permeability  of  2.75  ±  1.76  x 
10^  cm2  and  hydraulic  conductivity  of  15.1  ±  9.7 
cm/s,  values  analogous  to  a  gravel  deposit.  The 
relatively  high  coefficient  of  variation  (69%)  was 
attributed  to  large  variation  in  the  structure  of  the 
field  deposits  tested  as  indicated  by  borehole  re¬ 
sistance  logs. 

Ice  piece  size  and  distribution 

Little  information  regarding  ice  piece  size  and 
distribution  in  jams  is  available  in  the  literature.  A 
search  of  the  CRREL  Ice  Jam  Data  Base  (IJDB)  in 
June  1999  identified  some  information  on  ice  piece 
size  in  less  than  0.2%  of  the  records.  It  is  clear  from 
the  available  descriptions  that  damaging  ice  jams 
can  occur  when  parent  ice  thickness  ranges  from 
very  thin  (about  6  in.)  up  to  very  thick  (4  ft).  Brief 
estimates  of  floe  size  and  distribution  were  given 
for  two  breakup  jams  and  one  freezeup  jam.  The 
only  known  quantitative  size  distribution  is  for 
brash  ice  (Fig.  15).  The  data  were  developed  from 
a  photograph  taken  on  the  St.  Clair  River  in  Feb- 


Effective  Diameter  (cm) 

Figure  15.  Distribution  of  ice  piece  sizes  observed  on  the  south  channel  of  the  St.  Clair 
River,  25  February  1987.  (From  Daly  and  Arcone  1987.)  Ice  piece  diameter  is  ex¬ 
pressed  here  as  equivalent  area  circles  for  measured  ice  pieces . 


23 


ruary  1987  (Daly  and  Arcone  1989).  The  wide  vari¬ 
ability  in  piece  size  given  in  this  figure  is  also  evi¬ 
dent  in  qualitative  descriptions  contained  in  the 
IJDB  of  ice  piece  size  provided  for  breakup  jams. 
For  example.  White  (1996):  "The  ice  pieces  near 
the  toe  appear  to  be  thermally  grown  ice,  four  to 
eight  inches  thick,  with  occasional  thicker  pieces. 
The  pieces  are  neither  exceptionally  large  or  small, 
and  are  in  the  range  of  two  to  six  feet  in  diameter, 
with  the  mode  closer  to  the  smaller  size."  White 
(1993)  and  Tuthill  (1993)  also  made  observations 
of  piece  size  and  distribution  for  two  ice  jams  that 
formed  on  the  Black  River  in  Coventry,  Vermont, 
in  1993.  The  first  jam  formed  in  early  January  and 
was  made  of  blocks  generally  8  inches  to  10  inches 
thick,  ranging  up  to  18  inches.  This  jam  froze  in 
place  and  was  followed  by  a  later  jam  in  March, 
for  which  ice  pieces  averaged  1  to  3  feet  by  0.5  ft 
thick.  Lever  et  al.  (1997)  reported  that  ice  piece  size 
near  the  toe  of  breakup  jams  ranged  from  three  to 
four  times  the  ice  thickness,  similar  to  estimates 
by  Wuebben  and  Stewart  (1978). 


RELIABILITY  CONSIDERATIONS 


In  selecting  an  appropriate  range  of  values  for 
use  in  modeling,  it  is  valuable  to  have  some  idea 
of  the  reliability  of  the  properties  described  above. 
The  coefficient  of  variation  (COV)  can  be  obtained 
from  the  distribution  of  known  values: 


COV  = 


c(x) 

E(x) 


(30) 


where  a(x)  is  the  standard  deviation  of  the  sample 
and  E(x)  is  its  expected  value,  or  mean. 

Generally,  sample  sizes  of  the  ice  properties  dis¬ 
cussed  above  are  too  small  to  develop  this  infor¬ 
mation;  the  exception  would  probably  be  ice 
roughness.  However,  knowledge  of  the  COV  for 
various  soil  properties  will  provide  an  indication 
of  the  order  of  COV  for  similar  ice  properties.  Harr 
(1987)  summarizes  the  COVs  for  a  number  of  soil 
and  ice  properties  and  provides  references  as  well 
(see  Table  5).  He  notes  that  a  COV  less  than  about 
25%  indicates  a  level  of  reliability  acceptable  in 
most  engineering  applications,  while  a  COV 
greater  than  about  40%  indicates  unreliability. 

Using  the  information  in  this  table  as  a  rough 
guide,  we  could  consider  that  present  knowledge 
of  ice-influenced  velocity,  ice  roughness,  porosity 
of  frazil  ice  deposits  in  the  laboratory,  and  perhaps 
intrinsic  permeability  of  frazil  ice  in  the  labora¬ 


Table  5.  Representative  coefficients  of  variation 
(COV)  based  on  data  from  Harr  1987  and  some 
of  the  ice  properties  discussed  in  text. 


COV 

Property  Material  (%) 


Porosity  (p ) 

Soil 

10 

Specific  gravity  (y) 

Soil 

2 

Internal  angle  of  friction  (({)) 

Soil 

10 

Internal  angle  of  friction  for  sand  (<})) 

Soil 

7 

Internal  angle  of  friction  for  gravel(tj)) 

Soil 

12 

Cohesion  (c) 

Soil 

40 

Thickness 

Ice 

17 

Flexural  strength 

Ice 

20 

Crushing  strength 

Ice 

13 

Flow  velocity 

Ice 

33 

Ice  roughness  («t) 

Ice 

35 n 

Porosity  (frazil  ice  in  laboratory) 

Ice 

19b 

Intrinsic  permeability  of  frazil 

Ice 

32c 

(laboratory) 

Intrinsic  permeability  of  frazil 

Ice 

69d 

(field) 

Piece  size 

Ice 

100* 

Calculated  from  data  presented  in  White  and  Daly  (1997). 
Calculated  from  data  presented  in  Andersson  and  Daly  (1992). 
cWhite  (1991). 
dWhite  and  Lawson  (1992). 

Calculated  from  data  presented  in  Daly  and  Arcone  (1989). 


tory  is  reasonably  reliable.  It  is  clear  that  additional 
work  is  required  before  the  intrinsic  permeability 
of  frazil  ice  in  the  field  reaches  an  acceptable  level 
of  reliability.  The  same  is  true  for  the  coefficient  of 
internal  strength  (p)  and  field  values  of  porosity, 
whether  for  breakup  jams  or  freezeup  jams.  Based 
solely  on  experience  with  soils,  one  would  expect 
that  values  of  porosity  and  friction  angle  for  ice 
could  be  developed  with  some  reliability,  while  it 
might  be  difficult  to  achieve  the  same  level  of  reli¬ 
ability  for  ice  cohesion.  It  is  likely  that  further 
study  could  result  in  reasonably  reliable  estimates 
of  frazil  particle  size;  however,  it  is  unlikely  that 
the  same  can  be  developed  for  ice  piece  size. 


SUMMARY  AND  CONCLUSIONS 

An  overview  of  ice  jam  hydraulics  (juxtaposi¬ 
tion,  undercover  deposition,  and  the  formation  of 
an  ice  jam  from  shoving  and  internal  collapse) 
identified  a  group  of  variables  important  in  hy¬ 
draulic  and  physical  processes,  namely 

•  Ice  cover  roughness  (nj) 

•  Cohesion  (Q) 
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•  Coefficient  of  internal  strength  (|i) 

•  Angle  of  internal  friction  (<J>) 

•  Porosity  ( p ) 

•  Shear  force  on  the  underside  of  the  ice  cover 

M 

•  Frazil  particle  diameter  (d)  and  fall  velocity 
coefficient  (Q). 

The  ice-jam-related  literature  was  searched  for 
references  to  these  variables,  values  of  which  were 
collected  and  reported  in  the  present  study.  Based 
on  these  data,  recommended  values  are  shown  in 
Table  6.  At  the  present  time,  ice  cover  roughness, 
whether  for  frazil  ice  or  breakup  jams,  is  the  only 
variable  of  those  listed  above  that  is  known 
with  any  reliability.  Much  work  remains  to  be 
done  to  determine  reasonable  values  for  the  other 
variables,  especially  the  coefficient  of  internal 
strength.  Safety  issues  inherent  in  making  field  ob- 


Table  6.  Recommended  values  of  hydraulic  and 
physical  properties  affecting  ice  jams. 


Property  Range  of  values 


Ice  cover  roughness  («{) 

Freezeup  ice  cover 

Freezeup  ice  jam 

Breakup  (thin  ice) 

Breakup  (thick  ice) 

0.010-0.060 

0.020-0.10 

0.020-0.10 

0.035-0.15 

Cohesion  (Q) 

Freezeup  ice  jam 

Breakup  ice  jam 

960-1200  Pa 
0-100  Pa 

Coefficient  of  internal  strength  (p) 

1. 0-2.0 

Angle  of  internal  friction  (<J>) 

Frazil 

Rubble 

20°-45° 

40°-60° 

Porosity  (p) 

Freezeup  ice  jam 

Frazil  deposit 

Breakup  ice  jam 

0.35-0.45 

0.3-0.6 

0.35-0.8 

Shear  force  on  the  underside  of 
the  ice  cover  (t) 

Calculated  from 
other  data,  about 
10-20  Pa 

Frazil  particle  diameter  ( d ) 

1-10  mm 

Frazil  fall  velocity  coefficient  (Cf) 
5-mm-  to  1-cm-diameter  particles 

1.0 

servations  have  hindered  the  collection  of  field 
data  for  breakup  jams.  Advances  in  remote  mea¬ 
surement  technology  may  be  invaluable  in  this 
regard. 

Based  on  experience  with  soils,  one  would  ex¬ 
pect  that  values  of  porosity  and  internal  angle  of 
friction  for  ice  could  be  developed  with  some  reli¬ 
ability  in  the  future,  while  it  could  be  difficult  to 
achieve  the  same  level  of  reliability  for  ice  cohe¬ 
sion.  Reliable  estimates  of  frazil  particle  diameter 
and  fall  velocity  coefficient  could  be  achieved  with 
careful  study  The  primary  obstacle  in  developing 
such  estimates  appears  to  be  the  difficulty  in  de¬ 
vising  measurement  methods  for  a  material  close 
to  its  melting  point. 

The  literature  review  also  provided  information 
on  other  hydraulic  and  physical  properties  of  ice 
jams  that  are  of  use  to  researchers  and  engineers 
interested  in  ice  jam  hydraulics.  These  are 

•  Anchor  ice  growth 

•  Jam  erosion 

•  Ice  cover  velocity 

•  Coefficient  of  ice  loss 

•  Thermal  effects  on  ice  jam 

•  Permeability 

•  Ice  piece  size  and  distribution. 

Recent  work  in  the  area  of  anchor  ice  growth 
indicates  that  a  Froude  number  greater  than  0.2, 
and  preferably  larger  than  about  0.40,  is  required 
for  anchor  ice  growth.  Permeability  has  also  been 
the  subject  of  recent  attention,  particularly  in  the 
case  of  grounded  jams,  where  porous  flow  may 
be  important  in  determining  water  levels.  More 
work  needs  to  be  done  in  the  area  of  jam  erosion 
and  thermal  effects  on  jams,  particularly  for 
freezeup  jams.  Additional  data  on  the  coefficient 
of  ice  loss  could  be  quite  useful  in  determining 
jam  volume. 
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